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A  nonlinear  parametric  adjustment  model  with  n  observations  and  u  parameters 
can  be  interpreted  geometrically  as  the  description  of  a  u-dlmenslonal  "model 
surface"  embedded  in  an  n  dimensional  flat  "observational  .space".  Two  geomi'tric 
concepts  have  boon  Idontifjcd  for  the  least-sqnaros  resolution  of  this  model  it'" 
first  concept  is  linked  to  the  0au.'“;,  fopin  ».f  ‘Lj  ouifaco.  which  is  tlie  most  direct 
'  ,  ji;  li.o  parametric;  model;  however,  the  solution  is  Iterative  in 

nature.  The  second  t;c)nc:ppt  is  linked  to  the  func.-tlonal  form  of  the:  surface,  and 
it  leads  to  a  non  iterative  solution. 
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The  Gauss  form  of  t  iu>  surface  is  pi-eseiifetl  as 

X  ■-  X  ( u  )  ,  r  1  , 2 . 11  .  rt  1  . 2 . 11  . 

r  (X 

where  x  are  the  space  coorciiriates  and  u  are  the  surface  coordinates  of  (loints 
lying  in  the  surface.  Upon  Identifying  the  observations  with  the  observational 
space  coordinates  of  the  point  Q,  the  least-squares  solution  is  siM'n  to  correspond 
to  the  point  P  in  the  model  surface  such  that  tlie  distance  QP  is  a  minimum. 
Accordingly,  the  vector  PQ  is  orthogonal  to  the  model  surface. 


The  Ga.uss  form  can  be  transformed  into  the  functional  form  of  the  surface: 

F^(x‘  )  -  0  ,  r  =  1,2 . n  .  L  -  1 . 2 ,  .  .  .  . n- u  , 

where  the  superscript  L  distinguishes  ii  u  functions  "F"  one  from  anotlier.  Tliis 
formulation,  where  n-u  represents  the  number  of  degrees  of  freedom,  can  be  linked 
to  the  condition  method  of  adjustment  with  a  very  special  set  of  functions  "F". 

A  readily  perceived  advantage  of  the  functional  form  is  that  the  n  n  gradients  of 
these  "F"  are  vectors  orthogonal  to  the  model  surface,  which  is  instrumental  in 
expressing  the  final  vector  PQ. 

Both  concepts  make  use  of  the  Taylor-series  expansion  containing  the  first- 
and  higher-order  partial  derivatives  of  the  observables  with  respect  to  the 
parameters.  In  the  Gauss -form  approach,  the  derivatives  can  be  used  at  most  up  to 
the  third  order.  In  principle,  such  a  truncation  of  the  Taylor  series  results  in 
the  necessity  to  Iterate  the  solution.  In  the  functional -form  approach,  the  number 
of  useful  terms  in  the  Taylor  series  has  nc  limitation,  and  can  be  chosen  sucli  that 
the  final  outcome  Is  determined  with  sufficient  accuracy  in  a  single  step.  This 
approach  is  the  subject  of  a  detailed  study. 


The  role  of  eventual  nonlinear  constraints  in  both  geometrical  approaches  is 
to  restrict  the  model  surface  to  a  "model  subsurface",  but  otherwise  no  conceptual 
modification  takes  place.  Accordingly,  the  resulting  nonlinear  adjustment  model 
contains  a  reduced  number  of  parameters.  From  this  point  on,  the  treatment  of  tlie 
new  model  is  Indistinguishable  from  the  treatment  of  th('  original  model  without 
constraints , 
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1.  INTRODUCTION 


A  series  of  studies  linking  an  adjustment  model  to  a  geometrical  setup  with 
tensor  structure  and  notation  began  with  fltlaha,  19S4].  However,  that  report 
treated  a  linear  or  llneai'ized  adjustment  model;  a  nonlinear  model  would  In 
general  enta’1  an  Iterative  process.  A  more  recent  study,  [Blaha,  1989aJ  , 
contained  a  geometrical  development  of  an  accelerateo  itei'atlve  process  taking 
advantage  of  the  second-  and  third-order  partial  derivatives  (at  the  most)  of 
the  observables  wltti  respect  to  the  parameters,  in  addition  to  the  usual  first-^ 
order  partial  derivatives  featured  by  the  linearized  model.  The  present  study, 
by  contrast,  describes  a  non-iterative  approach,  where  advantage  can  further  be 
taken  of  the  fourth-  and  possibly  higher-order  partial  derivatives.  In  addition 
to  its  theoretical  interest,  this  approach  enables  one  to  resolve  problems  where 
fourth-order  partial  derivatives  are  significant  and  cause  difficulties  in 
alternate  adjustment  schemes. 

The  approach  is  based  on  the  Tay lor-series  expansion  of  the  observables  in 
terms  of  the  parameters,  and  on  the  geometrical  concept,  generalized  to  higher 
dimensions,  that  the  gradient  of  the  potential  at  a  given  point  is  orthogonal  to 
an  equipotential  surface  at  tliat  point.  We  emphasize  that  whenever  an 
adjustment  is  mentioned,  it  is  the  least -squares  adjustment  that  we  have  in 
mind.  In  geometrical  terms,  this  implies  the  generalized  notion  of  the  shortest 
distance  from  a  surface,  and  of  the  associated  orthogonality.  Indeed,  the 
least-squares  principle  has  provided  the  motivation  for  interpreting  adjustment 
models,  linear  or  nonlinear,  with  or  without  constraints,  in  a  generalized 
context . 

The  mathematical  background  needed  in  relating  nonlinear  parametric 
adjustment  to  geometry  is  presented  in  Chapter  2.  Parts  of  this  material  are 
common  to  some  of  the  previous  reports,  and  have  appeared  in  Chapters  1  and  2  of 
[Blaha,  1989a],  and  in  Chapter  2  of  [Blaha,  1989b],  However,  their 
recapitulation  is  useful  for  the  development  of  a  new  approach,  referred  to  as 
the  functional- form  approach,  in  the  body  of  the  present  study.  The  "functional 
form"  derives  its  name  from  tht;  manner  in  which  we  express  the  "model  surface"; 
the  latter  is  In  general  a  hypersurfaco  representing  the  adjustment  model  in  its 
geometric  analogy.  The  resulting  algorithm  has  been  implemented  on  a  personal 
com()uter.  and  verified  using  the  adjustment  problem  of  Appendix  A  in  [Blaha, 
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1989a],  formulated  in  terms  of  a  tKlrd-det^ree  polytiomial  iu  fout  variables.  The 
outcome  of  this  verification  is  presentftd  in  Apfiendix  A  herc'in. 

The  method  developed  for  non-iterative  adjustment  of  a  nonlinear  parametric 
model  can  be  used  whetht'r  or  not  there  exist  constraints  anionf;  the  pai-amet  ers . 

In  the  adjustment  context,  the  presence  of  s  constraints  can  be  used  to 
eliminate  s  parameters.  In  the  geometrical  context,  the  s  constraints  restrict 
tile  u  dimensional  model  surface  to  an  embedded  us  dimensional  "model 
sut.'surface  "  .  With  the  exception  of  the  model  subsurface  takiiig  place  of  the 
original  surface,  the  geometrical  concept  remains  ui'-'Hered.  Indeed,  from  this 
point  on.  the-  tr(;atment  of  the  new  model  is  indistinguishable  from  tlu'  treatment 
of  the  original  model,  whether  by  the  method  developed  herein  or  by  the  methods 
of  [illaha,  1989a].  The  rt'striction  of  the  model  surface  to  tlu;  mod(’l 
subsuiface,  and  Its  algorithmic  realization,  are  the  subject  of  Ap[)(Midix  B, 
which  represents  an  expanded  version  of  Sections  3.1  and  3.2  in  [Blaha,  1989b]. 


2.  ELEMENTS  OF  THE  MATHEMATICAL  BACKGROUND 


T!u'  pai’iinii't  r  i  f  .id  j  ns  t  mtMi  I  modi'l  imcIi  nf  t  lii‘  olisi'rvahlcs  in  li'rms 

of  iiiiinmc  t  crs  ,  whoi'o  the  sfruolnrt*  linkintf  •  li<‘  two  pioups  of  v.ii- i  .ili  I  cs  is,  in 
general,  nonlinear.  Tlie  number  of  observables  is  denoted  by  n  arui  the  number  of 
parameters  by  u,  where  n  must  be  greater  tlian  u  for  an  adjustment  to  take  place. 
In  matrix  notation,  t.he  adjustment  model  is  written  as 

I."'  -  Fix'*)  . 

wtiere  i.^  and  x'*  are  the  sets  (column  vectors)  of  adjusted  observations  and 
adjusted  parameters,  respectively.  Sucli  estimates  are  always  considered  to  be 
consistent  with  the  least -scpiares  principle. 

In  a  standard  adjustment  approach,  a  nonlinear  adjustment  model  is  subject 
to  the  Taylor-series  exfiansion  based  on  an  initial  set  of  parametric  values,  x'*. 
The  terms  in  the  second  and  higher  powers  of  the  parametric  corrections  are 
neglected,  resulting  in  the  well-known  (linearized)  observation  equations.  In 
matrix  notation,  the  latter  are  expressed  by 

V  =  A  X  +  L  , 

where  A  is  the  design  matrix,  X  =  X'*-X**  is  the  column-vector  of  parametric 
corrections,  V-L^-L^  is  the  column- vector  of  residuals,  and  L  =  L°-L**  is  tlie 
CO  1  umn- vector  of  constant  terms,  wltli  I,**^F(X**)  representing  observables 
consistent  wltli  the  initial  set  of  parameters,  and  1,^’  containing  the  actual 
observations.  The  linearized  model  is  subjected  to  the  least  squares  criterion 
T 

V  P  V  =  minimum  . 

where  F  is  the  weight  matrix  of  observations.  This  criterion  leads  to  the 
formation  of  the  familiar  normal  equations. 

If  the  original  adjustment  model  is  nonlinear,  the  resolution  of  the 

linearized  model  does  not  yield  the  final  answers.  The  process  is  usually 

repeated  with  new,  updated  parameters  and  the  corresponding  changes  in  A  and  L. 

However,  the  variance  lovariance  matrix  of  observations,  L,  as  well  as  the 

weight  matrix  P.  idofited  as  P  L  V  arc  constant  Thus,  the  matrix  of  noimal 
T 

equations,  N-A  I'A,  cliang.es  only  due  to  A.  and  the  column  vector-  icpiesent  i  ng 

T 

the  right  hand  side  of  normal  equations,  U^  A  PI.,  changes  only  due  to  A  and  1 
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The  edmim  I  a  t.  i  oti  of  tlie  uixiated  p.i'ametric  vaUn's  tlii’oinjli  a  ni-'.v  X  tequirt's  llie 
‘ermation  <tiul  the  i  livers  i  an  (if  a  new  N  in  eaeii  iteration,  oi  a  ma  iiemal  i  ea  1  1  y 
equiv’aJent  procedure.  Wtien  X  becomes  sufficiently  close  to  zero  the  iterative 
process  is  terminated.  As  its  by-product,  the  latest  matrix  N  '  is  adopted  as 
the  variance  covarlaiK'e  matrix  of  adjusted  pa’’amelers. 

The  functional  relationship  between  the  observable.^  and  the  parametc'rs 
lends  itself  to  a  geometrical  interpretation  and  treatment  Involving  spaces  and 
surfaces  generalized  to  higher  dimensions.  In  particular,  the  formulation 

( u* )  ,  -  2  2 . n  ,  a  1,2 . u  , 


represent  1  iig  tile  paramet>jc  adjustment  model,  can  be  llnktul  to  t  tie  (iauss  form  of 
a  surface  in  relation  to  the  surrounding  space,  wlutrt'  x'  are  t  h('  spare 
coordir.at  ‘s  and  u*  are  the  surface  coordinates.  The  Gauss  form  of  a  two- 
dimension-il  surface  (u--2)  embedded  in  a  three  dimensional  flat  space  (n^3)  is 
desrribeu,  togc'ther  with  two  other  forms,  in  Chapter  6  of  [llotine.  1969].  In 
[Blaha,  1934],  lioth  the  n-dimeiis  lonal  "observational'  space  and  the  u- 
dimensional  ’mode]’’  surface  were  considered  flat.  The  litter  was  thus  in 
reality  a  tiyperpJane.  Although  the  model  surface  is  now  intrinsically  a  curved 
spa(;f.‘,  the  surrounding  space  is  again  flat,  and,  as  is  shown  below,  its 
coordinate  system  is  characterized  by  a  constant  metric  tensor. 

p 

In  den.iting  the  n  observables  by  x  ,  r'1.2 . n,  and  the  u  unknown 

p.iram(He-s  by  u^ .  a^l.2 . u,  wo  can  represent  a  nonlinear  parametric 

adjustment  model  by 


[■ 

X 


r ,  a 
X  ( u  ) 


All**  +  (1/2)  n\  Au*  Au^ 

o  a  a$ 

*  (i  '6)  <Ii'  Au*  Au^  Au^  1  ... 


tia) 


Au 


a 


a 


u 


a 

u 

o 


(1b) 


"  r  r  Of 

wiierr.  represents  an  initial  set  of  parameters  and  x^=x  (u^)  represents 

I  lie  I  bservalile.s  consistent  with  this  set.  The  lower-case  Roman  indices  range 

froni  I  to  n,  and  flie  lower  case  Greek  indices  range  from  1  to  u.  Tensor 

symijoiisni  implies  the  summation  convention  over  the  dummy  (repeating)  indices 


In  tiie  geometrical  context,  the  first  ecpiality  in  (la)  represents  tlie  Gauss 
form  of  a  u  dimensional  surface  embedded  in  an  n  dimensional  space.  The  surface 
i I'lidowi'd  with  ttie  luiordinate  system  (u**),  a- 1,2 . u  and  is  referred  to  as 
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the  model  sitrface,  and  the  space  ji  endowed  with  the  .-nordiiiale  system  lx'}, 
r=!,2,.,.,  ,  and  is  ref(?ri'ed  to  as  tiie  oliservat  Iona  1  space.  ihe  second  eijuality 


in  (1,  '  is  the  Taylof  series  expansion  of  x  f i om  ttie  "initial"  point  P  lyinf.  in 


the  model  surface,  whose  model  -surface  coorditiates  are  u  and  whose 

o 

observational  -  space  coordiiiates  are  The  notation  identifying  liie  partial 

derivati''es  at  P,  such  as  dx*/du**=A*  .  3 ‘'x  *  /  d  ii**  i)  Ml*  ,  etc..  Is  s-lf 

a  a(S 

evident.  The  actual  observations  can  be  thought  of  as  describing  the  point  ' 


the  oliservat  i  ona  1  space,  which,  due  to  measurJrg  errors,  does  not  lie  in  ttu! 
known  model  surface.  The  task  at  hand  consists  in  determining,  from  tlie 


ohsi  tx'ed  point  Q,  a  model  surface  point  satisfying  tlie  1  east -squares  ra  Iterion 


In  the  adjustment  context,  the  variance-covariance  and  the  weight  matrices 
of  observations  depend  on  the  gu  ilty  of  measurements.  They  arc  indeper.dent  of 
ttie  adjustment  model,  of  the  initial  set  of  parameters,  of  the  outcome  of 
observations,  etc.  Ttvs ,  for  a  givcui  observational  nioiie  they  ;ir('  (constant  ,  In 
the  "traditional"  identification  of  fbiaha,  1984],  variance  covariance  matrices 
correspond  to  assficiated  metric  tensors,  and  weight  matrices  correspond  tc 
metric  tensors.  Acc'ird  1  ngl y ,  we  represent  the  varlance-covaria.acc  matrix  of 

p  g 

observations  by  tbe  oliservat iona 1 -space  associated  metric  tensor  g  ,  ami  the 
weight  matrix  of  observations  by  the  observational -space  metric  tensor 
state  that  both  tensors  are  Independent  of  the  form  of  the  model  surface,  oi  the 
initial  point  P.  of  the  ohserv(?d  point  Q.  etc.,  leading  to  the  simplification 


g^^  =  constant  (2) 

rs 

One  could  also  attribute  the  tensors  g  and  g^^^  to  the  point  Q  and  siate  that 
the  geometrical  setup  must  account  for  Q  located  anywhere  in  the  observational 
space.  In  turn,  (2)  implies  that  tlie  observational  space  mu.sl  he  flat. 


•  r 

If  the  set  (5x  denotes  the  coordinate  differences  between  the  observed 
point  Q  and  the  desired  model -surface  point  denoted  P.  It  corresponds  to  the 
neg.itive  residuals,  and  ttn'  least  sipiares  criterion  corresfioiids  to 


Ttie  quadratic  form  (3)  represents  the  square  of  the  liistance  between  Q  and  P 
Therefore,  the  desired  " 1  east  -  squares "  point  P  must  be  the  fovit  point  of  the 
straight  line  dropped  orthogonally  from  Q  onto  the  model  surface.  We  note  that 
if  any  other  adjustment  condition  wei e  used  in  lieu  of  the  1  east  -  squares 


cf  i  lor  i t'li ,  llie  minimum- lUst  aiice  property  (3)  would  not  exist  and  the  geometric 
tensoi  ial  treatment  of  th('  ad  justment  thfory  would  probably  be  much  mor(!  complex 
i  f  I'.o I  i  mposs  idle. 

A  conve:; : 'Mil  approach  fur  resolving  nonlinear  least  squares  problems 

;:onsists  in  using  an  Isomorpiiic  geometrical  setup  with  tensoi'  sti'ucture  and 

notation.  Such  a  link  is  highlighted  by  the  consideration  that  the  least 

squares  criterion  glve-s  rise  to  a  minimum  distance  property  Amoiig  tie  basic 

1  orr espondences .  the  number  of  observations,  n.  and  the  number  of  parameters,  u, 

define  thf  dimensionality  of  the  observational  space  and  of  the  model  surface, 

respectively  Since  the  constant  variance  covariance  matrix  of  observations,  Z. 

rs 

(orrespiinds  to  the  associated  metric  tensor  g  .  and  the  weight  matrix  of 
observations,  adopted  as  Z  corresponds  to  the  metric  tensor  g  .  the 
riiisei  .Mt  ictoil  space  is  endowed  with  a  coordinate  system  {x  }  sucl;  that 

i-  constant  .  g' ^  -  constant 

s  r 

i  he  sf't  [/’  (  f  ai  tual  observations  corresponds  to  the  set  x^  observat i onal - 
space  coordinate-  descrlbl.ig  the  point  Q.  All  possible  sets  of  adjusted 
ohserval i ons  (subject  to  no  criterion)  correspond  to  the  Gauss  form  of  the  model 
surface  endowed  with  a  Coordinate  system  {u*}.  which  has  already  appeared  above; 

X  =  x(u),  r-1,2 . n,  a=].2 . u. 

T  e  final  set  of  adjusted  parameters.  corresponds  to  a  particular  set 

ii'^  nf  Ri.,  it  1  sur^'ace  cooruinates  describing  the  least-square;s  point  P.  The  set 

III  iniMal  parameters,  ,  corresponds  to  the  set  u*  of  model -surf  ace 

cd'w  d  in<iles  describing  tlie  initial  point  P.  The  final  set  of  parametric 

(  !i:  I  Cl  I  i  (IMS  ,  X,  then  corresponus  to  Au**^u*-u*;  these  quantities  are  assumed 

to  Ilf  sni<i  1  1  .  The  final  set  of  adjusted  observations,  L^"F(X^).  corresponds  to  a 
r  r  cx 

particular  set  x  =x  (u  )  of  otiservational  space  coordinates  describing  the 

least  squares  point  P.  The  initial  point  P  is  described  by  these  coordinates 

as  .x'  .x'(u*),  reflecting  its  counterpart  L*^=F(X°).  The  set  of  negative 

o  o  j  o  r  r 

f  oust, ant  terms,  -  L  =  i,  F(X  ),  corresponds  to  the  contravaria.it  vector  6x  =Xq 

v'  ,  wliile  tlie  sot  of  negative  residuals.  V-h^  I,^,  corresponds  to  the 

f  I 

'lint  Mvatiaiit  v  ctor  f5x'~-.x'  x'.  The  initial  di.-sign  matrix.  A,  which  in 

Q 

siaiidard  observation  equations,  V  AX«h.  relates  the  parametric  coriections  to 

r  r  (1 

Iff  lesiduals,  corresponds  o  the  design  tensor  A^  =  0x  /au  evaluated  at  P. 
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On  the  other  hand,  the  standard  adjustment  approach  does  not  have  equivalents  of 

the  above  0*  and  <I>*  .  which  form  thiee-  and  four  dimensional  arrays, 

a(3  a$r 

respectively ,  and  contain  the  second  and  the  third  order  partial  derivatives  of 
with  rr^spect  to  u*,  evaluated  at  P. 

We  further  comment  that  the  surface  associated  metric  tensor  at  P 
corresponds  to  the  variance-covariance  matrix  of  adjusted  parameters,  and  the 
"necessary"  vrrrslon  of  the  space  associated  metric  tensor  at  P  corresponds  to 
the  singular  variance-covariance  matrix  of  adjusted  observations.  Such 
correspondences,  as  well  as  the  definition  of  "necessary"  tensors,  were 
presented  In  Sections  2.2  and  4.2  of  [Blaha,  1984].  (A  better  term  might  be 
"restricted",  in  the  sense  of  a  restriction  to  the  hyperplane  tangent  to  the 
model  surface  at  P.  )  In  the  same  vein,  the  surfacu'  metric  tensor-  at  P 
corresponds  to  the  weight  matrix  of  adjusted  parameters,  and  the  "necessary" 
version  of  the  space  metric  tensor  at  P  corresponds  to  the  singular  weight 
matrix  of  adjusted  observations. 

These  notions  confirm  the  variance-covariance  propagation  law  and  introduce 
a  "weight  propagation  law",  and  they  can  be  extended  to  functions  of  adjusted 
parameters  as  well  as  to  functions  of  adjusted  observations.  In  [Blaha,  1984], 
however,  they  were  demonstrated  only  for  linear  models.  In  nonlinear  models  the 
above  tensors,  expressible  by  means  of  orthonormal  vectors  tangent  to  the  model 
surface  at  P.  could  be  said  to  describe  tangential,  or  first-order,  properties 
of  least-squares  estimates.  In  the  present  study,  the  derivation  of  the  metric 
tensors  corresponding  to  the  weight  matrices  will  be  by-passed,  and  only  the 
derivation  of  the  associated  metric  tensors  corresponding  to  the  desired 
varlanc(?-covarlance  matrices  will  be  carried  to  term.  This  also  implies  that 
the  formation  of  the  weight  matrix  of  observations  (by  Inverting  their  variance- 
covariance  matrix)  will  be  by-passed. 


3.  FUNCTIONAL-FOPM  APPROACH  TO  NONLINEAR  PARAMETRIC  LEAST-SQUARES  ADJUSTMENT 


3 , 1  Two  Geometrical  Concepts  Reflecting  the  Least-Squares  MetLodoloKV 

In  the  preceding  chaptei',  we  have  seen  that  a  nonlinear  parametric 
adjustment  model  with  n  observations  and  u  parameters  can  be  Interpreted 
geometrically  as  the  description  of  a  u-dlmens lonal  "model  surface"  embedded  in 
an  n  (iimensional  flat  "observational  space".  The  most  direct  interpretation  is 
linked  to  the  Gauss  form  of  the  surface: 

=  x*^(u*)  ,  r  =  1,2 . n  .  -  1,2 . u  , 

r  Of 

wlieie  X  are  the  space  coordinates  and  u  are  the  surface  coordinates  of  points 
lying  in  tin?  surface.  Upon  identifying  t  l»e  ol)servatlons  with  the  observational 
space  coordinates  of  the  point  Q,  the  least  snuares  solution  corresponds  to  the 
point  P  in  the  model  surface  such  that  the  distance  QP  is  a  minimum. 

Accordingly,  the  vector  PQ  is  orthogonal  to  the  model  surface. 

In  a  separate  representation  as  treated  In  the  present  study,  one  can 
transform  the  Gauss  form  into  the  functional  form: 

F^(x^)  =  0  ,  r  =  1 ,2 . n  .  L  -  1  ,2 . n-u  , 

where  the  superscript  L  distinguishes  n-u  functions  "F"  one  from  another.  This 
formulation,  where  n-u  represents  the  number  of  degrees  of  freedom,  can  be 
linked  to  the  conrlltion  method  of  adjustment  with  a  very  special  set  of 
(Tinctioii.s  (■'*'.  A  readily  pcrtcived  advantage:  of  th(‘  functional  form  Is  that  tbc 
n-u  gradient.s  of  F*’  are  vector.s  oil  hogoJia  1  to  the  model  surface,  which  is 
i  us  f  r  umeiit  al  In  expressing  the  final  vector  PQ. 

Both  of  the  above  concepts  can  be  applied  to  the  nonlinear  parametric 

least  squares  adjustment,  where  use  is  made  of  the  Taylor  series  expansion 

containing  the  first-  and  higher-order  partial  derivatives  of  the  observables 
with  res|)iMt  to  the  paiami'ters .  In  tin'  Gauss  form  approach  developed  In  [Blaha, 
lUliPaJ,  the  partial  deiivatlve.s  can  h*'  used  al  most  up  to  the  third  order.  In 

principle,  such  a  truncation  of  the  Taylor  series  results  in  the  necessity  to 

iterate  the  solution,  but  this  process  converges  much  faster  than  would  be  thf* 
case  with  the  linearized  model.  On  the  other  hand.  In  the  functional -form 
approach  developecJ  herein,  the  number  of  useful  terms  in  the  Taylor  series  has 

ft 


no  HmiUitioti  atul  can  be  cbosen  siicb  that  P  is  det erm i /led  with  sufficient 
accuracy  in  a  single  step. 

Independent  of  the  mode]  surface  formulation,  the  1  f'ast-squares  piinciple 
is  illustrated  in  fig.  i.  This  illustration  Is  limited  In  that  n-3  (the 
observational  space  is  three  dimension.il).  u  =  2  (the  model  surface'  is  twcj 
dimensional),  and  LI  (L=n-u).  The  unbarred  quantities  in  the  figure  are 
associated  with  the  known  "initial  point"  P.  while  the  baiTed  quantities  are 
associated  with  the  desired  "least-squares  point"  P:  both  P  and  P  lie  Ir.  th^ 
model  surface.  The  vector  u  of  the  orthonormal  triad  6,  j,  v  is  orthogonal  to 
t  h  i .'.urfacr'.  with  a  similar  construction  .applying  for  P  and  other  surface 
poi  lit  s  . 

In  the  absence  of  observational  errors,  the  "observational  point"  0  would 
lie  in  the  model  surface.  (In  the  algebraic  formulation,  the  system  of 
observation  equations  would  be  consistent  and  there  would  bo  no  need  for  an 
adjustment.)  From  the  geometric  point  of  view,  an  optimal  estimate  of  a 
hypothetical,  error-free  observational  point  is  that  point  in  the  model  surface 
which  is  the  closest  to  the  actual  point  Q.  Thus,  the  line  connecting  Q  with 
this  closest  point,  whlcfi  is  P  by  construction,  is  nec^essarily  orthogonal  to 
the  model  surface.  This  is  Illustrated  in  Fig.  1  by  the  collnearlty  of  the 
vectors  PQ  and  v. 

Since,  as  has  been  discussed  in  Chapter  2.  the  metric  tensor  of  the 
observational  space  Is  constant  due  to  the  nature  of  the  parent  adjustment  setup 
( t  hr  weigf't  matrix  of  observations  i.s  Independent  of  the  mathematical  model,  the 
initial  values  of  the  parameters,  the  numerical  outcome  of  the  observations, 
etc  ) ,  we  have 

g  =  g  =  constant  . 
sr  sr 

Accordingly,  the  space  Christoffel  symbols  arc  zero  and  the  space  covariant 
derivatives  are  simply  the  partial  derivatives.  The  above-mentioned  optimal 
estimate  has  been  seen  In  (3)  to  correspond  to 

<5x^g  3x’‘  =  minimum  . 

sr 

where  the  summation  convention  for  the  repeating  (dummy)  Indices  applies.  This 
relation,  associated  with  the  point  P.  expresses  tho  1  east -stpiares  criterion. 
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THREE-DIMENSIONAL  REPRESENTATION 


OF  ISOMORPHIC  LEAST-SQUARES  SOLUTION 


OBSERVATIONAL 

SPACE 


GIVEN  PROPERTY:  g  =  g  =  constant 

s  r  s  r 


LEAST-SQUARES  CRITERION  IN  TENSOR  NOTATION: 

”  s  ~  r 

<5x  g  Sx  =  ainlHUB 
sr 

VECTOR  (5x  IS  ORTHOGONAL  TO  THE  MODEL  SURFACE  AT 


DEPICTING  THE  LEAST-SQUARES  SOLUTION 


"31 


3 . 2  Overview  of  the  Functional  Form  Concept 


In  accordance  with  the  foregoing,  a  least-squares  prohiem  can  he  viewed  as 
the  task  of  drawing  a  geodesic  (a  straiglit  line  in  the  flat  observational  spare) 
through  a  known  observational  point  Q  orthogonally  to  a  given  model  siiifaiu'. 

The  intersection  of  the  geodesic  wltl  h(“  model  surface  marks  the  desired  least 
squares  point  P.  In  the  functional  form  approach,  this  task  is  solved  by 
designing  a  geodesic  orthogonal  to  the  model  surface  tlirougli  the  initial  point 
by  implicitly  evaluating  the  miss  of  the  "target"  Q  by  this  line,  and  by 
dis(>lacing  the  point  P  In  the  model  surface  in  such  a  way  that  the  new  gi'odesii' 
(again  orthogonal  to  the  model  surfa<e)  passes  throu(;)i  the  target. 

The  functional -form  approach  can  be  described  as  consisting  of  four  steps. 
In  the  first  step,  the  equations  of  the  model  surface  are  transformed  from  the 
Gauss  form  to  the  functional  form  via  the  Monge  form,  leading  to  the  expressions 
for  the  (space)  gradient  vectors  associated  with  the  model  surface,  and  their 
derivatives.  In  this  process,  the  familiar  geometrical  concepts  are  extended  to 
a  u  dimensional  surface  embedded  in  an  n -dimensional  si)ace .  The  space 
coordinates  are  divid(!d  into  the  group  encompassing  the  first  u  coordinates, 
identified  by  the  lower-case  Greek  Indices,  and  into  the  group  encompassing  the 
last  n-u  coordinates,  identified  by  the  upper-case  Roman  indices.  Accordingly, 
the  Gauss  form  of  the  model  surface, 

X*  =  x'^(u^)  .  r  =  1.2 . n  .  j3  =  1.2 . u  .  (4) 

where  u^  represents  the  (original)  model -surface  coordinates,  is  written  as 

x“  =  x“(u^)  ,  (Ga) 

x^'  =  x^''(u^)  .  (5b) 


jO 

If  the  set  of  partial  derivatives  dx  /du*^  forms  a  regular  (i.e.,  square  and 
nonsingular)  matrix,  (x*)  will  be  adopted  as  a  n«-w  s«‘t  of  model  surface 
coordinates,  in  addition  to  repi’esent I ng  the  first  u  space  coordinates.  The 
Inverse  transformation  between  the  two  model -surface  coordinate  systems  is 
symbolized  by 


n 


(3 


(3,  a. 
u  (x  ) 


(Ga) 


The  substitution  of  (6a)  Into  (5b)  yields  the  Monge  form  of  the  model  surface: 


■;  1 


{6b) 


L 

X 


)  . 


wIhm'o  tlie  index  L  (or  anotlier  upper-case  Roman  index)  attributed  to  f 
distinguishes  the  n  u  functions  one  from  anotiier.  Finally,  the  functional  form 
of  tiie  model  surface  follows  from  (6b)  as 

x-f(x)=F(xJ=0.  (7) 


where  the  index  L  attributed  to  F  has  a  role  similar  to  that  just  mentioned.  In 
a  three-dimensional  space  containing  a  two-dimensional  surface.  (7)  would 
correspond  to  one  equation  F(xM=0.  such  as  F(x.y.z)=0  in  space  Cartesian 
coordinates,  whereas  in  the  general  case  being  developed  the  number  of  equations 
is  DU. 


In  general,  the  functional  form  based  on  (7)  is  represented  by 

F^'(x'  )  =  constant  ,  (8) 

where  "constant"  represents  a  set  of  n-u  constants.  The  zero  value  of  these 

constants  identifies  the  model  surface,  as  is  apparent  from  (7).  Here  again,  L 

(or  another  upper-case  Roman  index)  distinguishes  the  scalar  functions  of 

position,  N.  one  from  another,  so  that  its  role  is  to  account  for  all  of  the  n-u 

functions  N.  The  same  logic  is  utilized  when  describing  the  gradients  and 

higher-order  derivatives  of  these  functions.  Thus,  if  N  is  the  gradient  of  N, 

the  notation  symbolizes  the  entire  set  of  n-u  gradients.  The  advantage  of 

the  formulation  (8)  is  that  each  of  the  gradients  N  is  orthogonal  to  the 

hr  ^  r 

pertinent  surface,  in  the  sense  that  N^k  =0.  where  k  is  any  contravarlant 
vector  (i.e.,  a  vector  identified  by  Its  contravarlant  components)  tangent  to 
the  surface;.  This  orthogonality  property  is  crucial,  since,  as  we  have  seen, 
the  known  observational  point  Q  must  lie  on  a  line  drawn  from  the  least-squares 
point  P  orthogonally  to  the;  model  surface. 

In  addressing  this  task,  one  might  be  tempted  to  generate  a  family  of 
surface;s  in  space,  each  of  which  would  be  described  by  the  functional  form  (8) 
with  a  given  set  of  constants.  (When  these  constants  are  zero,  we  recover  the 
model  surface.)  We  could  stipulate  that  one  member  of  the  family  contains  the 
point  Q,  and  could  determine  its  set  N*'.  However,  this  avenue  presents  a  major 
difflfulty,  in  that  the  determination  of  intersections  of  straight  lines, 
orthogonal  to  the  model  surface,  with  the  new  surface  amounts  to  solving  a 
nonlinear  system  of  equations.  This  difficulty  lias  been  mentioned  merely  to 

1  ? 


illustrate  that  the  popular  concept  of  a  family  of  surfaces  does  not  advance  our 
task,  and  that  we  shall  have  to  resort  to  an  unorthodox  definition  of  a  surface 
as  Is  explained  below. 

The  secoi’d  step  of  the  new  approach  defines  a  special  "Q-surface" 
containing  th-;  point  Q.  However,  one  can  describe  this  surface  only  by 
construction  and  not  by  defining  equations.  In  particular,  a  point  of  the  Q- 
surface  is  generated  in  the  observational  space  from  a  chosen  model -surf ace 
point,  for  example  from  the  initial  point  P.  First,  we  conceive  of  n-u  straight 
lines  drawn  from  P.  which  are  orthogonal  to  the  model  surface  and  form  an  n-u 
dimension  il  iiy()iM  plane .  And  s(?cond,  this  hyperplane  is  Intersected  by  a 
straight  line  which  is  orthogonal  to  it  and  passes  through  the  point  Q.  The 
intersection,  denoted  P.  generates  a  point  of  the  Q-surface.  In  considering  the 
correspondence  between  a  general  point  P  in  the  model  surface  and  its 
counterpart  P  in  the  Q  surface,  we  confirm  that  the  poijit  Q  Itself  is  the 
counterpart  of  P, 

The  above  construction  can  be  11  Justrated  in  the  familiar  contoxl  of  a 
three-dimensional  space  containing  a  two-dimensional  surface,  where  only  one 
straight  line  perpendicular  to  the  model  surface  would  emanate  from  P.  This 
line,  representing  a  one-dimensional  "hyperplane",  would  be  Intersected  by  a 
straight  line  drawn  perpendicularly  to  it  from  Q,  which  would  generate  a  point 
of  the  Q-surface.  Such  a  situation  is  depicted  in  Fig.  2.  where  the  generated 
point  is  denoted  P  as  above,  and  v  is  a  vector  orthogonal  to  the  model  surface 
at  P. 


We  can  imagine  a  similar  situation,  illustrated  in  Fig.  3.  in  terms  of  a 
one-dimensional  "surface"  (l.e.,  a  line)  embedded  in  a  three-dimensional  space. 
There  would  now  be  two  vectors  emanating  from  P.  in  Fig.  3  denoted  v'  and  v  "  . 
which  would  be  perpendicular  to  this  "surface"  and  would  generate  a  "hyperplane’ 
(here  a  plane);  such  vectors  are  referred  to  as  vectors  v.  The  straight  line 
dropped  perpendicularly  from  Q  onto  this  "hyperplane"  would  then  generate  a 
point  of  the  Q-"surface",  denoted  again  P. 

The  construction  of  P  shows  that  the  vector  PP  is  the  orthogonal  projection 
of  the  vector  PQ  onto  the  n-u  dimensional  hyperplane  as  described  above,  and  it 
ensures  that  regardless  of  the  location  of  the  point  P  in  the  model  surface,  the 
Q-surface  passes  through  the  point  Q.  This  construction  has  the  advantage  of 


utilizing  a  set  of  gradients  (represented  by  the  vectors  v).  which  can  be 
easily  expressed  from  the  functional  form  of  the  model  surface  (as  obtained  In 
the  first  step),  and  of  projecting  Q  orthogonally  onto  the  hyperplane  generated 
by  these  gradients,  which  is  a  linear  operation. 

In  analogy  to  the  definition  of  the  (new)  model-surface  coordinates  x**,  in 
the  third  step  we  adopt  the  first  u  of  the  n  space  coordinates  x  as  the  Q- 
surface  coordinates,  denoted  ic*  for  easy  reference,  and  express  their 
variations  in  terms  of  variations  of  the  model -surface  coordinates  x**.  Since 
the  spatial  positions  of  points  in  the  Q-surface  are  linked  to  the  sets  of 
gradients  at  the  corresponding  points  in  the  model  surface,  the  variations 
in  x  entail  the  differentiation  of  with  respect  to  x“.  This  differentiation 
(a  part  of  the  first  step)  is  facilitated  by  the  relative  simplicity  of  the 
functional  form  (7).  The  mapping  of  the  variations  of  x®  into  the  variations  of 
x**  can  be  expressed  by  the  Taylor  series,  similar  to  a  coordinate  transformation 
in  one  and  the  same  surface. 

Such  a  series  can,  in  turn,  be  inverted  analytically,  and  the  inverted 
series  can  be  used  to  describe  the  variations  of  x^  in  terms  of  the  variations 

^  (X 

of  X  .  The  important  characteristic  of  the  process  resulting  in  the  original 
and  the  inverted  series  is  that,  again,  only  linear  operations  are  involved. 
Since  both  the  point  P  (generated  from  the  initial  point  P)  and  the  point  Q  are 
known  in  all  n  coordinates,  and  since  both  lie  on  the  Q-surface.  their 
coordinate  differences  AX*  are  known.  In  using  the  inverted  series,  one  can 
then  determine  the  unknown  model -surface  coordinate  differences  Ax*  between  the 
points  P  and  P. 

After  the  determination  of  the  coordinates  x*.  the  remaining  coordinates, 
x^,  follow  from  (6b).  Thus,  insofar  as  the  observational-space  coordinates  are 
concerned,  the  geometrical  problem  of  determining  P  is  solved.  In  the 
adjustment  context,  this  represents  the  solution  of  adjusted  observations  but 
not  adjusted  parameters.  As  the  fourth  and  last  step,  the  related  geometrical 
task,  l.e.,  the  determination  of  the  original  model-surface  coordinates  of  P, 
can  be  accomplished  via  the  initial  inverse  equation  (6a).  This  relation,  like 
others  encountered  above,  is  expressible  by  the  Taylor  series  expansion  from  the 
initial  point  P,  and  can  thus  be  resolved  by  means  of  linear  operations. 
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OF  A  ONE-DIMENSIONAL  Q-SURFACE 
4E-DIMENSIONAL  MODEL  SURFACE 


4.  DEVELOPMENT  OF  THE  FUNCTIONAL- FORM  APPROACH  IN  FOUR  STEPS 


4  .  1  FiiSt-Sten  Df.rKation 

r  r  6 

Iti  (4).  Ihi;  model  s'.irfaco  j^iven  in  the  Gauss  form  via  x  =x  (u  ).  This 

form  is  known,  as  are  the  partial  derivatives  listed  here  to  the  third  order: 


3x  /3u  =  A 

a 


^2  r  a  ^ 
d  X  /  3  u  3  u 


il  „  ,  3  X  13 n  3u  3u 

a/S 


®  .  ( 9a , b , c ) 

afir 


Each  of  these  sets  is  partitioned  into  two  sul)sets: 
,  r  , .  o  .  „r  ,  „p  jc  .  ,„r 


« 


{A^.  A*^}  . 
a  a 


if  s  (4)^  (J)*^  1  . 

a^r  a$r‘  a^r 


( 10a . b , c ) 


where,  in  agreement  with  the  previous  convention,  the  lower-case  Greek  indices 
range  from  1  to  u,  the  upper-case  Roman  indices  range  from  u+1  to  n.  or, 
equivalently,  from  1  to  n-u,  and  the  lower-case  Roman  Indices  range  from  1  to  n . 

In  consulting  (6b)-(8),  we  define  (space)  scalar  functions  of  position 
as 

-tL  L  0!.  L  L,  (X .  K 

N  -X  -f(x)sx  -x(u(x)). 

and  stipulate  that  along  the  model  surface  they  take  on  the  constant  value  of 

P 

zero.  Their  gradients  in  the  coordinate  system  (x  >  follow  as 

N^'  s  3N^/ax‘'  H  {aN^/3x^,  =  {-3x^/3x'',  ish  . 

r  R 

which  we  write  in  the  form 


N^  =  aN^/3x^  =  -ax^73x^ 
p 

N„  H  3N  /ax  =  6  : 

K  K 


(lia) 

(lib) 


we  keep  in  mind  that  the  quantities  x  in  (11a)  and  the  equation  above  it 
represent  f^.  functions  of  the  first  u  space  coordinates  as  introduced  in  (6b). 
The  nonzero  second-  and  third-order  partial  derivatives  are 


H  3^N^/3x^3x^  =  -0  x^/3x^3x^  , 
PU 


(12a) 


.,L  ^3.,L.-  P,  r 

N  £  3  N  /3x  3x  3x 
,J(IT 


3^x“/ 3x^dx^3x^ 


1  7 


(12b) 


If  we  adopt  {x**}  as  a  new  set  of  model  surface  coordinates,  the  Gai’^^s  form 
of  tlip  surface  can  be  written  as  x^=h^(x*),  representing  {x^,  x^)  ^  {x^.  f*^(x*) } . 
The  new  i)arttaJ  derivatives  (listed  to  the  third  order)  become 

ax‘'/ax“  =  (ax^/6x*,  ax^’/ax®}  =  {6^.  .  (i3a) 

a  a 

a^x''/ax“ax'®  =  id^  x^/ dx°‘ dx^ ,  a^x^/ax“dx^)  =  {o,  ,  (i3b) 

«p 

a^ x^/ ax® ax^ ax^  =  {a^x^/ax®ax^ax^ ,  0^x^/ax®ax^3x^}  =  {o.  }  .  (i3c) 

apr 

It  should  be  mentioned  that  these  (as  well  as  the  previous)  partial  derivatives 
are  consideied  at  the  point  P  of  the  model  surface.  This  understanding  applies 
also  for  the  subsequent  paragraphs. 


Next,  we  differentiate  (15)  further: 

a^x''/ax''ax^  5  3/ax^(ax^/ax'")  =  a/ax^[(ax''/au^)(au^/3x*)] 


The  right-hand  side  of  this  equation  Is  composed  of  two  terms,  the  first  of 
which  is  obtained  with  the  aid  of  9  /  3x  ^  ( 3  3ii^e  ii^) '  3  u  )  . 

Thus ,  we  have 

3^x'',  3x^9x^  =  (  3  ^x*  /  3  u^3  u  '  )  (  3u^/ 3x^ )  (  3  u^Z  3  ) 

^  (9:;''z9u^)(9^u^Z3x'^9x^)  .  (17) 


According  to  (13b),  the  left-hand  side  above  is  according  to  (9b) 

together  with  (10b),  and  (14a),  the  factors  in  tlie  first  term  on  the  right-hand 

side  are  0^  =(0^  ,  },  R^,  and  R^;  and  according  to  (9a)  together  with  (10a), 

p7  pr  P7  k  r 

and  (14b),  the  factors  in  the  second  term  on  the  right-hand  side  are 
r  H  S 

A„  =  (A'„,An)  and  A  Equation  (17)  is  seen  to  represent  two  equaliti<-s. 

p  p  P  ^  ^ 

Upon  the  contraction  with  R  and  the  use  of  (16a),  the  first  equality  yields 

P 

A'^  -•  R'^0^  R^R^  .  (18a) 

KT  H  PJ  K  Z 


With  this  result  as  well  as  (,.ub),  the  second  equalit;,  yields 


N 


L 

KT 


(N^'nJ 

P  $7 


fi7  K  t 


(18b) 


The  changes  in  i.idices  (two  in  dufmy  indices,  one  in  a  free  index)  which  have 
entered  these  two  formulas  need  no  elaboration. 


Similar  to  the  differentiation  of  (15'  which  has  lead  to  (17)  and  finally 
to  (18a. b),  we  now  differentiate  (17)  with  respect  to  x^.  Straightforward 
algebraic  manipulations,  where  advantage  is  taken  of  the  fact  that  th“  partial 
derivatives  commute,  again  result  in  two  equalities.  The  first  reads 


KTW 


Rh< 

^  pjo  K  X  W 


o'"  (R^A^ 

K  xw 


R^^ 

r  Q)K 


)) 

(D  KT 


(19a) 


and  the  second  reads 


=  -  [(N^D^  +  )  R^r’^K*^ 

KZii)  p  pj6  pyo  K  X  ut 

<■  (t/'n^  )(R^A^  +  R^a’'  +  R^a’'  )] 

p  fir  (9r  K  xm  X  UK  w  KT 


(19u) 


,L 


In  passing  from  (19a)  to  (19b).  R  has  been  replaced  by  N 'fl^  *0.  . 

'  p  P7  A  u  P  P7  P7 

the  same  as  when  passing  from  (18a)  to  (18b).  and  has  been  replaced 

.  .,L  .p  .1, 

by  N  i|i„  c  ■ 

ii  prs  PrS 


In  summary,  the  first-  .  second  .  and  tliiid  orZ  t  partial  derivativ  -s  ol  t  tie 
original  model -surface  coordinates  u*  with  r«?spect  to  the  new  model  surface 


coord  i  ria  t<‘s  as  outlined  in  (14a. b.c)  are  t’iven  by  ( 3  f>a or  (J6a'',  (38a),  and 

(19a),  rcsfioc  t  i  VO  1  y .  And  tin-  first-,  second-,  and  third  order  partial 
deri'.'allves  of  the  scalar  functions  N*'  with  rc-spect  to  the  new  model -surf ace 
coordinates  x**,  initially  presented  In  (11a)  and  (12a, h),  are  givcii  by  ^16b), 
(18b),  and  (19b),  respectively. 


4  2  Q  Surface  Determination 

The  d(  t  ei  ni  i  Hit  i  on  o)  t  tie  Q  surfacr;  will  now  bi  described  in  deta.i.  Since 
any  ctui.sen  point  of  ttie  model  surface  generates  a  unique  corresponding  point  in 
th  ■  Q  siiid  .K  e,  w-  ne<?d  only  to  express  a  typical  vector  <5x ,  such  as  Illustrated 
in  Fig.  2  and  3.  The  coot  lavar iant  and  cc.ariant  components  of  this  vector  are 

p 

denoted  6>:  and  respectively,  linked  together  in  a  familiar  relation: 

6x'  --  g''^6x  .  '20) 

s 

wtiere  is  the  associated  metric  tensor  (here  constant).  The  contravar iant 

p 

com[)onents  (5x  of  the  vector  <Sx  (illustrated  in  Fig.  1-3)  represent  the  basic 
known  quantities.  The  n-u  gradients  nK  in  Fig.  2  and  3  represented  by  the 
vectors  v,  comprise  anotfier  set  of  known  quantities  associated  with  P.  Since 
Sx  is  orthogonal  to  the  mode)  surface.  It  must  be  expressible  as  a  linear 
coml/ i fi.i 1 1  on  (jf  ttie  vectors  v.  In  covariant  components,  we  tints  have 

Sx  -  t/’c„  .  (21) 

where  r  are  n-u  coef f i c;ients :  the  index  M  (or  another  upper  case  Roman  index) 

M 

attt  Ibnted  to  c  acts  only  in  a  '  ^ountable"  role,  similar  to  that  explained  below 

ecjii.itlon  (8)  in  conjunction  witli  the  functions  N.  Thus,  c^^  is  not  a  covariant 

tensor  itnt  is  contracted  with  in  accordance  wi»h  the  summation  convention. 

s 

On  ttie  otiier  hand,  dx'  is  stipulated  to  be  a  vector  tangent  to  the  model 
suiftKc  at  f’ ;  it  is  thus  oi’ttiogonal  to  the  v<*ctors  v.  This  vector  is 
constrained  by  the  requirement  that  when  it  is  added  vectorially  to  <5x,  the 
result  is  dx  (see  also  Fig.  2  and  3).  In  contravariant  components  (and 
similarly  in  rovarian  components),  one  then  has  tlie  tensor  equation  at  P: 

<5x  ' '  <5x^  -  (S,x ' 


2  0 


(22) 


r  f' 

where  i5x'  and  6x  are  unknown  so  fai-.  The  above  tangential  property  is 
reflected  In 

N^<5x''’  -  (  dx^'  '  <5)1^ )  0  ,  (2.-1) 

r  t' 

where  Hie  "O'*  on  the  rig^ht  hcind  side  represents  ii  -  u  zero  elements.  At  this 

p 

stage'  w(?  use  (20)  and  (21)  In  expressing  <5x  in  terms  of  the  unknown 
coefficients  and  tlie  known  gradients,  so  that  (23)  becomes 

N^g''®N'^c„  -  N^'6x"  .  (24) 

r  s  M  r 

If  the  gradients  of  N  for  each  L  (or  M)  are  considered  as  column-vectors  of 

r  s 

n  elements  grouped  into  the  matrix  R*  of  dimensions  nx(ri-u/,  the  components  g 
are  grouped  Into  the  pos  1 1 1  ve- def  1  n  1  te  (symmetric)  matf'ix  g  of  dimensions  nxn. 
the  components  are  grouptul  into  the  column  vector  Ax  of  dimensions  n  1  .  and 

the  coefficients  syml)olized  by  c  are  grouped  into  the  col umii -vector  c  of 

M 

dimensions  (n-u)><l.  equation  (24)  rends 


(R*^gR*)c  =  R»^Ax  . 


(24'  ) 


The  transformation  from  (24)  to  (24')  does  not  reveal  the  usual  correspondence 

between  tensor  contractions  and  matrix  multiplications  as  we  would  expect  it  if 

L  and  M  were  "regular"  tensor  Indices,  since  In  such  cases  we  would  treat  the 

contravar iant  (upper)  index  as  the  first,  and  the  covariant  (lower)  index  as  the 

second.  However,  in  their  "countable”  role.  L  and  M  are  written  in  the  upper 

position  in  N^.  etc.,  for  convenience,  and  could  as  well  have  been  utilized  in 
r 

the  form  ^  ,  etc.  In  fact,  in  view  of  matrix  transcriptions  (which  are 
completely  unnecessary),  such  an  alternative  form  corresponds  to  the  above 
construction  of  the  matilx  R* .  and  n-sults  in  (24')  by  virtue  ol  the  summation 
convention  used  in  (24). 

bM 

Next,  an  array  quantity  B  is  introduced,  containing  (n-u)x(n-u)  scalar 
invariants  associated  with  P.  according  to  the  definition 

-  N^'g^^N^  .  (25) 

r  s 


In  matrix  form,  this  can  be  transcribed  as 


B  =  R’»^gR* 


(25'  ) 


whore  B  is  a  posi live- def Inl te  (symmetric)  matrix  of  dimensions  (n-u)x(n-u). 

duo  to  the  faf:t  that  is  pos  i  t  i  ve- dof  ini  te  (symmetric)  and  that  K*  has  the  full 

column  rank  n-u  (it  contains  n  ii  independent  Kfadlent  vectors  generating  an  n-u 

dimensional  hyperplane).  W(;  now  consider  an  array  quantity  Q  defined  through 

KL 

the  Kronecker  delta  such  that 

b“  =  ,  (26) 

which,  in  matrix  notation,  would  correspond  to 

QB  =  I  ,  Q  -  b'^  --  (K*^gR*)'*  .  (26') 


In  contracting  (24)  with  Q  and  recalling  (25)  and  (26),  we  obtain 

KI« 


where  =c  has  been  utilized.  In  matrix  notation  this  would  read 

KM  K 

c  =  QR*^<5x  =  (R*^  gR*)'^R*^(5x  . 


(27'  ) 


Finally,  with  the  aid  of  (21)  and  (27),  the  contravariant  vector  (5x ^  from 
equation  (20)  Is  expressed  as 

.,r  rs.,K  rs-.K,.,  .,L .  8)  , 

«  -  6  -  g  •  I 

which,  in  matrix  notation,  would  be 


<5x  =  gR*c  =  gR*  (R*^gR*)'^R'»^(5x  , 


(28-  ) 


where  the  n  components  of  <5x  have  been  grouped  into  the  column  vector  5x . 
As  a  matter  of  interest,  we  remark  that  due  to  the  general  relation 


1  u 

k  =  g  k_ 


k  =  g  k 
s  sr 


r  s  L 

the  set  of  (covariant)  gradient  vectors  could  have  been  transformed  by  g 

into  a  set  of  contravariant  vectors,  which  could  have  then  been  grouped  into  the 

matrix  R  of  dimensions  nx(n  u).  This  amounts  to  writing 


R  =  gR*  ,  R*  =  g*R  , 

upon  the  realization  that  g"(g*)  where  g*  groups  the  components  of  ggj.- 
Accordingly,  eejuations  (27')  and  (28')  v^ouid  have  been  formulated  in  terms  of  R: 


c 


T  -IT 
(R  g*R)  R  g*Sx 


Sx  =  R  (R^g*R)  ^  R^g*<5x  , 

which  have  a  form  akin  to  that  for  the  parametric  corrections  and  the 

observational  corrections,  respectively,  in  the  parametric  least-squares 

adjustment,  with  R  corresponding  to  the  design  matrix,  g*  corresponding  to  the 

weight  matrix,  and  <5x  corresponding  to  the  vector  of  (reduced)  observations. 

However,  the  original  formulation  transcribed  as  (27')  and  (28')  is  far  more 

advantageous  in  the  present  context,  since  the  initially  known  quantities  are 

the  gradient  sets  forming  R*  and  the  observational  variance-covariances 

r 

forming  g. 


Clearly,  the  matrix  transcriptions  and  the  adjustment  analogies  have  been 
incidental  to  the  present  development,  and  could  have  been  left  out  altogether. 
However,  their  presentation  may  be  useful  for  their  Instructive  value, 
especially  in  considering  that  the  adjustment  topics  have  been  traditionally 
exposed  in  matrix  notation.  In  closing,  we  recapitulate  the  outcome  related  to 
the  Q-surface  determination  in  a  concise-manner : 

<5x''  =  .  (29a) 

where 

c„  =  Q„,  nJ;<5x"'  .  (29b) 

K  KL  m 

and  where 


„LM  „L  rs  M 
B  =  N  g  N 
r  s 


(29c) 


Third  and  Fourth  Steps 


In  its  geometric  analogy,  nonlinear  parametric  adjustment  has  beeti  seen  to 

correspond  to  the  task  of  drawing  a  geodesic  in  an  n-dimensional  observational 

space  orthogonally  to  a  u  d i mens  1 ona 1  model  surface.  Tlie  relative  simplicity  of 

this  ttisk  stems  from  the  fundamental  property  that  the;  observational -space 

metric  tensor  Is  constant.  Accordingly,  the  observat lonal -space  coordinate 

r 

system,  symbolized  by  (x  ).  can  be  thought  of  as  oblique  Cartesian,  and  a 

r 

coordinate  set  x  of  a  givtm  point  can  be  though  of  as  a  contravariant  position- 
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p 

vector  p  .  This  Interpretation  allows  us  to  write  what  may  appear  as  a  mixture 

r  r 

of  vector  components  (here  <5x  )  and  coordinates  (here  x  ); 

<5x''  =  x^(Q)  -  x'^(P)  .  (30) 

where  the  letters  in  parentheses  denote  the  points  in  question.  Before 
commenting  further,  we  introduce  the  notational  conventions  as  follows. 

(i)  The  point  Q,  corresponding  to  the  outcome  of  observacions ,  is 
considered  fixed,  and  its  coordinate  sets  are  henceforth  attributed  the 
subscript  Q  instead  of  being  identified  ty  "(Q)". 

(ii)  The  point  P  is  considered  to  be  a  variable  point  of  the  model  surface. 

and  its  coordinate  sets  are  henceforth  attributed  no  identification;  thus, 

r  V 

inst(!ad  of  writing  x  (P),  we  will  write  simply  x  . 


(iii)  The  point  P  of  ttie  Q-surface.  corresponding  to  the  point  P  of  the 
model  surface,  could  have  been  identified  by  the  coordinate  set  written  as 

r  *  p 

X  (P),  similar  to  the  original  notation  x  (P)  seen  in  (30).  However,  due  to  the 

letters  In  parentheses  being  dropped,  we  replace  "x"  by  "x"  for  any  point  of 

the  O-surface;  thus,  x*^(P)  and  x^(P)  are  replaced  by  x''  and  s'',  respectively. 

This  convention  applies  also  for  subsets  of  the  first  u  or  the  last  n-u 

ct  oc  Ij  L 

coordinates,  such  as  in  x  and  x  ,  or  x^  and  x^.  When  referring  to  the  fixed 
point  Q  (rather  than  to  the  variable  point  P  or  its  counterpart  P),  we  use  the 
subscript  Q  as  stipulated  in  item  (i),  i.e.,  we  write  x  and.  eventually,  x 

r  V  V 

or  xj. 

With  the  above  cor, '.nt  ions ,  equation  (30)  is  written  as 
r  r  r 

dx  =  Xq  -  x  .  (31 ) 

^  r  r 

If  and  X  are  interpreted  as  contravariant  position  vectors  belonging  to  P, 
then  (31)  represents  a  tensor  equation  at  P.  In  fact,  due  to  the  character  of 

P  p 

(x  )  alluded  to  above  (oblique  Cartesian  coordinates),  the  coordinate  set  x 
can  be  Identified  with  the  contravariant  position-vector  associated  with  P. 
Although  X  is  identified  with  the  contravariant  position  vector  associated 

P 

with  Q.  the  character  of  (x  >  allows  this  position  vector  to  be  freely  parallel 
transported  to  P  or  any  otlier  location;  the  "free"  parallel  transport  Implies 
that  no  change  in  the  components  has  taken  place.  Similarly,  we  have 


..r  .r  r 
(5x  =  X  -  X 
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which  can  again  be  interpreted  as  a  tensor  equation  at  P,  The  geometrical 
objects  representing  the  contravariant  vectors  and  <Sx*^  are  6x  and  <5x, 
referred  to  simple'  as  vectors,  which  are  illustrated  in  Fig.  2  and  3.  In  the 
following,  we  will  consider  only  a  subset  of  (32)  encompassing  the  first  ii  out 
of  n  comn<ir>e»ts  tHo  rpinfinp  ttamely 


..vtt  .a  a 
ox  =  X  -  X 


(32'  ) 


is  not  a  tensor  equation,  but  merely  an  equality  for  the  selected  components. 
This,  however,  does  not  detract  from  its  usefulness  in  the  present  development. 

We  use  (32  '  )  in  I  h('  f orm 
a  .»a 

St  =  X  +  6x  .  (33) 

from  which  we  can  deduce  that  x^  is  a  function  of  x”.  In  fact,  from  (28)  or 

(X  r*  rs  ots 

(29a),  6x  follows  from  <5x  upon  replacing  g  by  g  : 

6x“  =  ■  <34) 

PS  ots 

First  of  all,  g  is  constant;  this  property  affects  g  above,  as  well  as^^Q 

LM  L 

via  B  (see  equation  29c).  From  (lla.b)  it  is  apparent  that  is  a  function 

r 

of  X  .  Finally,  fix  is  written  as 

fix'"  2  (fix^,  fix^)  ,  (35) 


where,  due  to  (31),  one  has 


fix^  =  x{f  -  x^ 


.  M  .M  M 

fix  =  Xq  -  X 


(35'  ) 


M  M  u 

Here  x^  and  x^  are  constant,  while  x  is  a  function  of  x  as  was  already  stated 

in  (6b),  in  conjunction  with  the  Monge  form  of  the  model  surface.  It  then 

follows  that  fix®,  and  therefore  x®,  are  functions  of  x**  (see  equations  34 

and  33).  A  similar  statement  can  also  be  made  with  respect  to  the  "full" 

version  of  (33).  corresponding  to  equation  (32). 

We  are  now  in  a  position  to  present  the  concept  of  the  determination  of  the 
least  .squar(!S  point  P  based  on  tiin  role  played  by  the  Q  surface.  This  point 
is  regarded  as  the  location  of  the  variable  point  P  which  entails  x”'*!!®, 
l.e.,  for  which  the  Q-surface  points  P  and  Q  coincide. 


(a)  The  crucial  element  of  the  least-squares  resolution  is  the  set  of  u 
coordinates  x**  of  the  variable  point  P  when  the  latter  coincides  with  P. 
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(The  romainiiiK  n- u  space  coordinates.  x*‘,  can  lx;  found  via  the  Monge  form  of  the 
model  surface  as  presented  In  equation  6b,  and  the  u  model -surface  coordinates 
u**  for  this  point  can  be  found  via  equation  6a.) 

(b)  The  set  x*  associated  with  the  variable  point  P  can  be  obtained  from  a 
cpt-  of  !•  variables  if  the  functional  relationship  between  the  two  sets  is  known. 
(If  this  relationship  is  formulated  through  sets  of  partial  derivatives,  the  set 
x**  can  be  expressed  by  the  Taylor  series  expansion:  such  a  formulation  is 
sufficient  to  consider  the  relatlonslilp  as  known.) 

(c)  The  just-mentioned  u  variables  are  represented  here  by  the  coordinate 

ck  r  r 

set  X  .  To  elaborate,  we  first  state  that  the  sets  x  and  <5x  are  known 

ot  r  r 

functions  of  x  .  Accordingly,  due  to  x  =x  e^x  .  we  can  also  state  that  the 

r  ot 

set  X  and  any  of  its  subsets  are  known  functions  of  x  .  A  simple  and  natural 

manner,  then,  in  whlcli  to  choose  a  suitable  set  of  u  variables  consists  in 

adopting  the  first  u  out  of  n  relations  represented  by  this  tensor  equation, 

namely  x^=  x**+ (5x** ,  which  has  appeared  as  (33). 

(d)  By  virtue  of  the  known  functional  relationship  between  the  sets  x**  and 
x**  associated  with  P  and  P.  respectively,  we  can  express  the  coordinate 
differences  Ax®  In  terms  of  Ax®  and  vice-versa.  Since  this  known  relationship 
reflects  the  fact  that  the  vector  PP  is  orthogonal  to  the  model  surface,  any  new 
pair  of  points  P  and  P  (found  upon  applying  the  coordinate  differences  to  a 
previous  nalr)  will  preserve  this  quality. 


(o)  We  next  consider  on<;  position  of  P  as  known,  but  keep  the  notation 

unaffected  by  this  designation.  Such  an  "initial"  point  is  given  by  the  model- 

Ot  r 

surface  coordinates  u  .  from  which  the  observational -space  coordinates  x  follow 

by  (4).  The  sets  dx^ .  Ax®,  and  X®  associated  with  this  point  are  then  known 

from  (31),  (34),  and  (33).  respectively.  Since  x®  is  known,  so  is  now 

Ax®Hx®-j(®.  This  leads  to  Ax®  and  thus  to  the  mode) -surf ace  coordinates  of  a 

new  point  P,  whose  corresponding  point  P  Is  Q  (due  to  the  definition  of  Ax®). 

Accordingly,  the  new  ;;oint  P  is  the  desired  least-squares  point  P. 


(f)  Since  the  (new)  model -surface  coordinates  have  been  defined  as  {x®}, 
and  similarly  the  Q-surface  coordinates  can  be  defined  as  (x®),  the  procedure 
described  in  items  (d)  and  (e)  transforms  coordinates  between  two  surfaces. 
This  approach  is  based  on  the  known  functional  relationship  between  the  two 
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kinds  of  coordinates,  dictated  by  the  known  correspondence  between  points  in  one 
and  the  other  surface  (points  "P"  and  points  "P").  We  comment  that  such  a 
procedure  is  akin  to  the  transformation  of  coordinates  in  one  and  the  same 
surface,  based  on  a  known  functional  relationship  between  the  coordinate  sets. 


The  implementation  of  the  above  concept  proceeds  along  the  following  lines. 

Since  the  observables  correspond  to  the  observat ional -space  coordinates  {x  )  and 

the  parameters  correspond  to  the  model -surface  coordinates  {u*}.  the  adjustment 

r  r  (x 

model  is  interpreted  in  the  geometrical  context  as  x  =x  (u  ).  appearing  in  (4). 

(X  r 

The  known  coordinates  u  of  the  initial  point  P  then  lead  to  x  .  The  first-, 

r  (X 

second-,  third-,  ...  order  partial  derivatives  of  x  (u  )  with  respect  to  the 

r  r 

model -surface  coordinates,  evaluated  at  P.  compose  the  arrays  A  .  0 


a&7 


The  design  matrix  A  of  the  standard  adjustment  setup  corresponds 


to  the  "design  tensor"  A  The  actual  observations  are  interpreted  to 
^  r 

correspond  to  x  .  and  the  variance-covariance  matrix  of  observations  is 

rs 

interpreted  to  correspond  to  g  ,  the  (space)  associated  metric  tensor;  the 

latter,  just  as  the  metric  tensor  g  .  is  stipulated  to  be  constant.  The  known 

rs 

^  r  r  r 

Xy  and  X  give  5k  by  (31),  which  then  leads  to  6k  according  to  (34). 

Since  a  close-form  relationship  between  x^  and  k°‘  via  (33)  and  (34)  would 
be  intractable  and  nearly  impossible  to  establish,  we  relate  the  two  sets 
through  their  partial  derivatives; 


*a  ,  «a  ,  ^  3  .a  .  _  jS 

A„  =  dx  /0x  =  <5„  +  0  5x  /0x  . 
p  P 

s  0  X  /0x  0x  =  0  <5x  /0x  0x  , 
07 


(36) 

(37) 


etc.,  which  are  based  on  (33).  These  arrays  can  be  expressed  from  (34).  where 

OtS 

all  the  quantities  on  the  right-hand  side  except  g  (constant)  are  functions  of 
x*.  to  be  evaluated  with  the  coordinates  of  P.  To  formulate  the  set  of 
coordinate  differences  Ax*  defined  in  item  (e)  above,  we  first  write 


dx'*^  =  <5x*^  -  dx^  =  , 

where  the  first  equality  is  (22),  and  the  second  equality  is  the  consequence  of 

p 

(31)  and  (32).  In  view  of  the  character  of  the  coordinate  system  {x  } ,  it  Is 

P 

not  surprising  that  dx'  ,  a  (space)  contravariant  vector  at  P.  is  expressible 
in  terms  of  coordinate  differences.  Upon  considering  the  first  u  out  of  n 
component. s  in  this  relation,  it  follows  that 
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(38) 


.  .a 

Ax  .  Xq 


a  .a  .»a 
X  =  ox  -  (5x 


We  remark  that  this  equation,  together  with  the  first  equation  of  (35'),  yields 

(38-  ) 


sax'*/  0x^ 


,  .a  fi , 

+  S6x  /3x  ) 
p 


which  Is  -A^  The  left-hand  side  of  (38')  is  one  of  several  expressions 
encountered  in  an  initial  formulation  of  0*  . 

|3r 


Equations  (36),  (37), 


would  allow  us  to  express  Ax'*  from  Ax'*  via  the 


Taylor  series.  However,  here  we  are  seeking  Ax'*  based  on  Ax'*.  The  reversed 
partial  derivatives,  needed  for  a  reversed  series,  are  denoted  as 


8x  /ax  = 

P 


-2  -« 

a  x  /ax  ax  =  a„ 


They  can  be  obtained  from  A* ,  0*  .  . 

0  07 

i n  R*^ .  A*  ,  ...  in  {16a),  (18a) . based  on  A* ,  0*  ,  .  .  . 

0  07  0  07 


upon  using  the  approach  that  resulted 

That  is  to  say. 


equations  (16a),  {18a),  ...  can  be  re-used  with  all  the  quantities  attributed 


the  symbol 


Thus,  the  array  R  is  determined  from 
P 


a'^R'^  =  <5® 

A  K  K 


(39) 


in  matrix  notation  this  implies  that 


[r;^]  =  ca“j'^ 


(39*) 


where  the  indices  are  now  meaningless.  Similarly,  the  array  is  given  by 


r^r"  . 

iSr  a  vr  0  7 


(40) 


The  above  equations  may  appear  simple,  but  they  are  greatly  complicated  by 

the  complexity  of  A® .  .  For  example,  we  can  show  that 

P  07 


-a  a  an  a  L 

“f!  ■ ^  ’’■0  '  “[.''e  ■ 


(411 


where 


e 


E  „  =  , 

00  K  u0 


G 


a 

L 

aK 


P'**^Q 


8 


KL 


0^  =  AK^ 


ar 


00 


aK 


l/'  =  =  N^Sx^  *■  Sx^  ; 

m  n 
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furthermore . 


B  .  or  [0^1(1  -  IB  1  . 

wlu'rc 

^KM  „K  rs.,M  ,.K  r<5  K<5  ,,K  (5M  KM 

B  =Ng  N  =  Ng  N-  +  c  Nj  +  N.g  +g 
r  s  7"  o  o  o 

In  which  we  have  used  (lib).  As  is  illustrated  above,  also  the  quantities  . 

.  ...  are  Involved  in  the  formation  of  A* ,  Q*  .  Fortunately,  these 

firi  /9  |9r  a  a 

quantities  have  been  found  in  conjunction  with  the  development  of  R  ,  A  ,  ... 

P  PT 

(see  equations  16b,  18b,  ...  ).  The  coordinate  differences  Ax  in  the  formula 
for  are  those  featured  in  (38). 

Tlie  derivation  leading  to  (41)  will  now  be  outlined.  We  begin  with  (34), 
which  is  transcribed  as 

Accordingly,  we  have 

a5i?"/ax^  =  (ap“'^/ax'®)Q^j^w^'  +  p“‘^[(aQ|^^/dx^) w^  +  (aw^/ax^) ]  , 

where 

r* 

-n<xK,-  /5  at.,K 

ap  /dxr  =  g  , 

aw^/ax^  =  N^^Ax^  +  N^(8ax^/ax^)  +  dSx^/dx^  =  N^.Ax^  : 

ufi  fl' 

the  last  outcome  is  due  to  tlie  cancellation  stemming  from  d6x^/ dx^  =  -6^  and 

aax^/ a x^  =  - ax^/ 0x^  =  .  The  expression  for  dQ  /dx^  follows  from  the 

P  MK 

differentiation  of  the  form  B  Q  =<5V.  where  this  equation  itself,  as  well  as 
MK  ^ 

the  symmetry  of  B  and  hence  of  Q  ,  are  also  utilized.  The  result  is 

K  L 

aq^j^/ax®  -  -(Gjq„  .  o“q^s)N®j  . 

Upon  using  the  previously  Introduced  notation,  we  obtain 

aax"/ax^  =  t“^e  ,  .  g“u^  , 

Pfi  L  (3 

and,  due  to  (36),  verify  (41). 


The  principles  that  have  led  to  (41)  can  be  further  expanded  to  yield 
0*  We  list  the  final  results  of  this  development  in  the  form  0^  6x^6x^ 
as  needed  for  the  contraction 


where 


A®  Ax^Ax’'  =  -R®0«  r'^r’'ax‘'^Ax^  =  -R®0^  dx'^dx^  , 
i9r  fi  $r  1^ 


,  $  ~B  _<o 
dx^  =  R^Ax  . 

(0 

In  a  form  suitable  for  programming,  where  vectors  such  as  dx®  first  contract 
higher-dimensional  arrays  in  order  to  lower  the  dimensionality,  we  have 


dx^ix^ 


^  P0  " 


Where 


p,  =  c  „dx®dx^  +  2d  dx®  , 
4>  K  <t>aB  K  <t)a 


“k  ■ 


.dx®dx^Ax^  -  dx®dx^  -  2N*^  dx®  ( d Ax^/ 3 x^ )  dx^  , 
ua$  a/9  fia 

where  0(5x^/8x^  is  given  explicitly  at  the  close  of  the  preceding  paragraph. 

The  desired  coordinate  differences  Ax®  are  obtained  with  the  aid  of  the 
Taylor  series: 

Ax®  =  R^Ax*®  +  (1/2)  A®  AX^Ax^  +  .  .  .  .  (43) 

P  pT 

It  is  Interesting  to  note  that  AX*  appears  not  only  in  the  direct  form  on  the 

right-hand  side  of  (43).  but  also  in  the  formation  of  R®.  A®  ,  ...  as  has  been 

P  P  7" 

illustrated  above.  This  stems  from  the  design  of  the  Q-surface,  involving 

implicitly  the  point  Q  already  in  the  construction  of  any  point  P;  the  more 

obvious  role  of  Ax®,  reflected  explicitly  in  (43),  is  that  of  "bridging  the  gap" 

between  P  and  Q.  Another  interesting  feature  of  the  series  (43)  is  that  when  It 

is  expressed  in  terms  of  the  initial  partial  derivatives,  the  linear  term 

involves  not  only  A®  but  also  fi®  ,  the  quadratic  term  involves  A®.  Q®  .  and 
a  a  P  a  ^  fir 

.  etc.  If  Ax  should  be  expressed  from  Au  .  the  linear  term  would  involve 

only  A®  (it  would  read  A®Au^).  the  quadratic  term  would  Involve  only  0®  .  etc. 
fi  fi  fir 

Since  the  coordinate  differences  Au®.  corresponding  to  the  parametric 

corrections,  are  important  quantities  to  be  resolved,  it  is  expedient  to  compute 

Fj  ot 

them  tiefore  completing  the  set  Ax  (by  its  subset  Ax").  From  the  subset  Ax 
('(jmpiited  via  (4.t).  we  obtain  Au®  as 


Au®  =  R®Ax^  +  (1/2)A®  Ax^Ax''  (1/6)0®  ,Ax'^Ax’'ax®  + 

p  pr  pro 
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whose  present  form  corresponds  to  that  of  (43)  for  the  reason  explained  prior  to 
(39).  We  can  express  (44)  alternately  in  terms  of  ,  0®  ,  ,  ...  as  is  seen 

0  0r  0r&  ^ 

below  in  (45a, b).  In  both  (44)  and  (45a. b)  the  linear  term  involves  only  . 

P 

the  quadratic  term  involves  only  A„  and  0.  .  the  cubic  term  involves  A®, 
a  «  0  0r  0 

,  and  <D„  etc.  To  present  all  three  series  to  within  the  third-order 

0r  0rS  ^  ^ 

initial  partial  derivatives  (i.e.,  <D„  _).  we  include  an  additional  term  in 

0r<5 

(44)  and  (45a)  as  compared  to  (43).  The  alternate  expression  (45a, b),  suitable 
for  programming,  is  written  as 


Au®  =  du“  -  (l/2)R®(n^j  du^di/)  -  (1/6)r“(®^  ^du'^du^du'^ 

fi  07  H  0r& 


3<  du^R^'o'^  du'^du^l 

07  n 


(4  5a  ) 


where 


.a  r>®  A 

du  -  R  Ax 
w 


(45b) 


Finally,  the  remaining  subset  to  be  resolved  is  found  from  the  simplest  of 
the  four  series: 

Ax^  =  a';Au^  +  (1/2)0^  Au'^At/  ^  (1/6)(I>^  .Au^Au’^Au^  +  ...  ,  (46) 

p  07  070 

which  corresponds  to  (la)  limited  to  the  last  n-u  coordinate  differences.  The 

model -surface  coordinates  of  the  least-squares  point  P  are  u°+Au**.  while  its 

r  r 

observational -space  coordinates  are  x  +Ax  .  where  the  former  corresponds  to  the 
adjusted  parameters,  and  the  latter,  to  the  adjusted  observations.  One  could 
verify  the  consistency  of  the  derivations  and  the  accuracy  of  the  computations 
by  re-evaluating  Ax*,  this  time  based  on  Au*: 

Ax*  =  A*Au^  +  (l/2)fl*  Au*^Au^  +  (1/6)  4)*  .Au’^Au^Au*^  +  ...  : 

P  p7  P70 

the  series  corresponds  to  (la)  limited  to  the  first  u  coordinate  differences. 

All  of  the  observational -space  coordinates  of  P  could  be  verified  by  using  (4). 

We  remark  that  the  parametric  resolution  via  (44)  can  be  by-passed,  and  the 
subset  Ax^  can  be  computed  by  means  other  than  (46).  Based  on  Ax*  found  in 
( 43 ) ,  we  have 


Ax**'  -  -N*^Ax*  (  1 /2 )  N*^  Ax*  Ax^  -  (1/6)  Ax*Ax^Ax^  - 

a  0(3  afir 


which,  with  the  notation  from  (45b),  is  developed  as 


Ax*^  =  A*^clu”  +  ( 1 /2 )  (N**' (0^  du^du^)  ♦  du^du^] 

a  '  /Sr 

(1/6)n’^[((Ii^  du^du^du**  -  sa'i  du^R^'o’^'  du'^dii^) 

/'  ^r<5  0r  ri  Kx 

+  <ti*^  du^du^du*^  -  3fl*^  du^R^CO*^  du*du^)]  +  ...  .  (46') 

0r<5  $r  r)  Kt 

The  expressions  written  here  in  parentheses  are  identical  to  tlieir  counterparts 

(X  K 

in  {45a):  their  storing  can  be  useful  when  both  Au  and  Ax  are  to  be  computed. 
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5.  VA!  lANCE  COVARIANCKS  AND  OTHER  ADJUSTMENT  QUANTITIES 


Vve  now  turn  to  tho  task  of  finding  the  associr.ted  metric  tensor  and 

the  "necessary"  (or  restrlctc'd)  associated  metric  tensors  ■it'd  H  ■  whet' 

tli(!  overbar  implies  tlie  evaluation  at  the  least  sipiaies  poin*  f’ .  The  first  of 
those  tensors  belongs  t  r)  the  u  dimensional  model  surface  embedded  in  the  ii- 
di Miens ’ ona 1  observational  space,  and  corresponds  to  the  var i ance -covar i ance 
matrix  of  adjusted  parameters.  The  second  tensor  be.ongs  to  the  observational 
space  restricted  to  the  u  limonsional  hyperplant'  tangent  to  the  model  surface  at 
P,  and  ('oi  resp  inds  to  t  htr  variance  covariance  matrix  of  adjusted  obsi'rvat  i  ons . 
And  the  third  tensor  belongs  to  tlu;  observational  space  restricted  to  the  n  u 
dimensional  subspace  orthogonal  to  this  hyperplane,  and  corresponds  to  the 
variance-covariance  matrix  of  residuals.  The  orthogonal  su! space  just  mentioned 
will  be  referred  to  as  an  n-u  dimensional  "complementary  hyperplane".  At  the 
point  P,  we  have 

-rs  -,rs  -„rs 
e  •-  B  +  g 

r  s 

where  g  is  the  observat ional -space  associated  metric  tensor  corresponding  to 
the  variance-covariance  matrix  of  observations.  Alt.iough  g' ^  is  considn  ed  to 

be  constant  for  any  location  in  the  observational  space  (see  Chapter  2).  the 

rs  “  rs  “ 

same  cannot  be  said  about  g'  and  g"  ,  which  depend  on  tl”'  location  of  P. 


The  present  approach  will  by-pass  the  formation  Oi  the  model -surface 

metric  tensor  at  P.  which  would  correspond  to  the  weight  matrix  of  adjusted 

parameters,  or.  equivalently,  to  the  matrix  of  nor-a]  equations  updated  for  P. 

Likewise,  It  will  by-pass  the  formation  of  g  =b  .  the  observat ional -space 

rs  rs 

metric  tensor,  which  would  correspond  to  the  (constant)  weight  matrix  of 
observations.  Due  to  the  fact  that 


i  .  g 

ap  a 


”S 


st 


which,  in  matrix  notation  are  equivalent  to 


ap' 


Ib 


rs 


fg 


St,  1 


it  follows  that  certain  mati ix  inversions  will  be  avoided.  For  the  sake  of 
simplified  notation,  the  desired  tensors  will  be  derived  for  the  initial  point 
and  the  resulting  formulas  will  then  be  applied  at  the  point  P.  We  thus 
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.  j  ,  .  ,.rs  , 

pi  0  0 «■  (Hi  lo  develop  a  .  g'  ,  and  g  .  where 


rs  ,rs  „rs 
K  H  +  g 


(47) 


rs  rs  (x3 

These  teuiors  will  be  derived  in  the  order  g"  ‘  .  g'  ,  and  a  . 


We  liegin  the  development  bv  reviewing  the  notion  of  the  "design  tensor" 
appearing,  e.g.,  in  Section  T  1  and  Appendix  B  of  fBlaha,  1989a].  This 
tensor  can  be  expressed  by  means  of  orthonormal  vec^’n'-s  £,  j,  ...  spanning  the 
1  porplane  tangent  to  the  model  surface  at  P: 

=  ax‘'/au®  =  »  j^j  t  . . .  .  (48) 

a  o  a 

r 

with  }  r-1,2 . n  repreo  ..nt  i  ng  the  observational -space  coordinates,  and 

(u*),  a-1.^  ...,u  representing  the  tangent -hyperplane  coordinates.  The 
contravar iant  components  of  a  vector  k  belonging  to  the  tangent  hyperplane 
transform  as 


r  r  cx 

k  =  A^k  .  (48') 

a 

The  design  tensor  (48)  can  serve  In  formulating  the  metric  tense  ^  the  tangent 
hyperplane,  * 


%  ■  *  b-'. 


aS  ,r 


(49) 


as  well  as  the  necessary  associated  metric  tensor  restricted  to  the  tangent 
hyperp 1 ane , 

r'o  r>o  T'  rtG 

(50) 


. ts  „r  s 
g  =  C  C 


.r  ,s 
J  J 


.r  ap  s 
A  d  A. 

«  B 


which  are  both  symmetric. 


If  the  i ndepe-ndent  variables  (u*)  represent  the  coordinates  for  both  the 

tangent  hyperplane  and  the  model  surface  its;'lf,  the  tensors  in  (48),  (49),  and 

150)  can  be  used  for  both  he^e  surfaces.  However,  such  equivalencies  hold  true 

only  at  P,  whose  model -surface  coordinates  as  well  as  tangent - hyperplane 

coordinates  are  u* ;  whereas  the  three  tensors  are  constant  in  the  tangent 
o 

hyinrplaiH,’ ,  their  values  in  the  model  surface  change,  in  general,  from  point  to 
p.iiit  (sec  Appendix  B  in  fBlaha.  1989a]).  Since  this  chapter  is  concerned  only 
witii  tip'  val'ies  at  one  point  the  conceptual  difference  between  the  model 
curtain  aiul  the  tangent  hyperplane  can  be  Ignored. 
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In  Chapters  3  and  4  we  have  seen  that  the  n-u  dimensional  complementary 


hyperplane  is  defined  by  n-u  gradients  (orthogonal  to  the  model  surface), 
L=1.2 . n-u.  Tlie  set  of  contravar lant  compoiients  of  these  n-u  gradients  are 

,Lr 


N 


rs  ^,1, 
g;  N 


(51  ) 


We  now  adopt  (x  ).  K=1.2 . n-u  as  the  coordinates  for  the  complementary 

p 

hyperplane,  and.  In  analogy  to  (48)  and  (49),  write  A  for  its  design  tensor 

K 

and  a  for  its  metric  tensor.  Thus,  in  analogy  to  (48'),  if  v  is  a  vector 
KL 

belonging  to  the  complementary  liyperplane,  we  have 
r  .  r  K 

V  =  A^v  . 


For  a  set  of  n-u  vectors  v,  whose  space  contravarinnt  components  are  grouped 
Li  r 

into  N  and  whose  complementary-hyperplane  contravarlant  components  are  grouped 
..LM 


into  N 


it  follows  that 


.Lr 


a"n^^ 


(52) 


R 


If  we  wish  that  for  each  vector  v  the  contravarlant  components  v  be  the  same  as 

r  r 

the  last  n-u  components  of  v  ,  the  last  n-u  coordinates  of  the  system  (x  )  must 

b  at  the  same  time  the  complementary-hyperplane  coordinates  (x  }.  (In  this 

case,  the  components  of  aJI  are  {A^,<5^).)  This  useful  provision  gives 

K  K  K 


s  (N^‘°,N^^} 


(53) 


Lr  Lr  LH 

where  N'  is  known  from  (51).  The  symbols  N  and  N'  do  not  represent  tensors, 

r  R 

but  only  a  collection  of  n-u  contravarlant  vectors  v  and  v  ,  respectively. 


We  next  define  the  symbol  such  that 


""""JK  ■  "kl" 


LM 

=  • 


or 


(54) 


The  combination  of  (52)  and  (54)  yields 
,,Lr 


Lr 


(55) 


(Upon  considering  -equations  53  ar.d  54,  this  confirms  the  parenthetical  statement 
from  the  close  of  the  preceding  paragraph.)  In  analogy  to  (49),  we  deduce 


ML 


I  c  Ip 
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where 


.,Js  ,,Ir  .,J  rs,,I  ,  ^IJ, 

N  =  N  N  =  N  g  N  =  B  (-0  . 

s  r  r  r  s 

hence 

“mi.  ' 

KL 

Since  the  associated  metric  tensor  a  is  such  that 

KL  ~K  ,  KL,  , 

a  ^ 

and  since 

>*  'JjK  -  'JkL®  '  \  "  ">KlI  ■  <®  '  ■ 

which  has  already  appeared  in  a  similar  form  below  (41),  we  deduce: 

a'<^  =  .  (58) 

In  fact,  by  contracting  this  formula  with  a  from  (56).  one  recovers  equation 

LN 

(57).  In  analogy  to  (50),  it  follows  that 
„rs  .r  KL.s 

g  =  Aj^a  : 

the  prime  has  been  replaced  by  the  double  prime  because  the  role  of  the  tangent 
hyperplane  has  been  taken  by  the  complementary  hyperplane.  Upon  using  equations 
(55)  and  (58)  in  conjunction  with  (54),  this  tensor  becomes 

g"^s  =  .  (59) 

PS 

We  can  take  advantage  of  the  quantities  already  found,  and  formulate  g" 
in  terms  of  Its  subsets: 


^.RS  _ 


From  ( 51 )  we  have 


^  SO 


the  quantity  has  appeared  below  (41)  directly,  while  the  result  for  N  has 

KM  PS 

participated  in  forming  B  (ibid.).  The  subsets  of  g"  are  thus  presented  as 


g  .  cjp"'*  l.g-‘"’l  . 


(60a) 


(60b) 
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(60c  ) 


„RS 


""'Ojn" 


NS 


<  ..SR, 

(=g  )  . 


From  (58)  we  observe  that 


.,RS 


RS 


(61  ) 


This  kind  of  equivalence  arises  whenever  a  hyperplane's  coordinates  coincide 
with  a  set  of  space  coordinates. 

rs  rs 

Having  found  g"  '  ,  we  obtain  g'  from  (47): 


g' 


rs 


rs  „rs 
g  -  g 


(62  ) 


rs 


For  the  first  subset  of  g'  .  It  follows  that 


_ 


,po 


=  -  o-p 


!  62  '  ) 


Next,  we  consider  (x  }  to  be  a  coordinate  system  for  the  tangent  hyperplane, 
identical  with  the  first  u  coordinates  of  the  system  {x  }.  In  denoting  the 
tangent-hyperplane  associated  metric  tensor  in  these  coordinates  by  in 
analogy  to  (61)  one  has 


(63) 


Whereas  3*^  Is  the  associated  metric  tensor  In  the  system  {x°'},  we  are 
Interested  In  a®^,  the  associated  metric  tensor  in  the  system  {u®}.  Due  to  the 
transformation  formula 

a  =  (0u  /ax  )S  (au  /ax  )  , 

where  au**/ax^  =  R®  has  already  been  employed  In  (14a),  equation  (63)  together 
with  (62 ' )  yields 


.  R«(gP"-  gJp'’''ir« 
p  No 


(64) 


rs  rs  ctS 

Thus,  at  P  the  desired  tensors  g"  ,  g'  .  and  a  are  given  respectively  by 
(59)  or  (60a-c),  by  (62),  and  by  (64). 

r» 

Next  we  consider  (51c‘  from  (28),  written  as 


<5x'' 


N  c 


K  ' 


(65) 


the  two  subsets  of  which  are 


6x^  =  P^^- 


K  ’ 


<5i? 


=  n"' 


(65'  ) 
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From  (65) .  the  covariant  components  are 


<5x  =  g  <5x^  =  N^r 

r  rs  r  l. 


where  the  last  expression  is  also  given  in  (21).  The  square  of  the  magnitude  of 
the  vector  <5x,  i.e.,  the  square  of  the  distance  between  P  and  P.  is 

-.2  o-r,.  v,Kr.,L 

<5s  =  <5x  <5x  =  N  N  c„c, 

r  r  K  L 

However,  since  according  to  the  equation  preceding  (56)  we  have  it 

r 

follows  that 


6s  =  c^B  Cj^  ^  c^B  . 

and  thus  (see  the  equation  below  57); 


6s^  =  c  W*^ 
K 


-  pg  -  pg  ~(xS 

The  desired  tensors  g"  ,  g'  .  and  a  are  obtained  from  their  unbarred 
counterparts  in  (59)  or  (60a-c),  in  (62),  and  in  (64),  respectively,  upon 
computing  all  the  pertinent  quantities  at  the  least-squares  point  P.  That  is, 
after  the  model-surface  coordinates  of  P  have  been  found  (denoted  here  as  u*) , 
the  algorithm  is  re-started,  but  the  computation  is  limited  to  the  quantities 

featured  below  equation  (41),  beginning  with  g“  The  computation  includes  also 

r  2 

the  contravarlant  vector  6x  in  (65)  or  (65'),  and  the  invariant  <Ss  in  (66). 

—  —r  “2 

At  the  point  P,  these  two  tensors  are  denoted  5x  and  6s  ,  respectively, 

representing  the  contravarlant  version  of  the  vector  ^  and  the  magnitude 

squared  of  this  vector.  In  the  adjustment  terminology,  6x^  corresponds  to  the 

-2 

minus  residuals.  -V.  and  6s  corresponds  to  the  least-squares  quadratic  form, 

V  PV.  All  of  the  above  tensors  a**^ .  dx*^ ,  and  6s^  are  verified 

numerically  in  Appendix  A. 


6.  DISCUSSION 


6 . ]  Relation  to  the  Condition  Method 

The  condition  method  features  adjustable  observables  without  any  adjustable 
parameters  present.  Its  geometrical  analogue  can  be  expressed  in  general  by 

F^(x^)  =  =  0  ,  (67a) 

where  has  been  Inserted  for  the  sake  of  comparison  with  the  outliru;  in 
Chapter  3.  In  fact,  from  equations  (7)  and  (8)  we  can  present  the  model 
developed  in  this  study  as 

H  -  f^(x“)  =  =  0  .  (67b) 

The  comparison  of  (67a)  and  (67b)  reveals  that  tlie  method  developed  herein  could 
be  regarded  as  a  special  case  of  the  condition  method,  in  the  sense  that  the  set 
of  functions  F^  has  a  very  special  form.  It  is  thus  more  conveniently  regarded 
(and  treated)  as  a  general  version  of  the  parametric  method. 

In  geometric  analogy,  the  parameters  have  been  seen  to  correspond  to  the 

model -surface  coordinates  {u®},  and  the  observables  have  been  seen  to  correspond 

r 

to  the  observational-space  coordinates  (x  ).  The  Gauss  form  of  the  model 
surface  represents  the  model  relationship  between  the  two  sets  at  a  general 
point ; 

r  r,  a, 

X  =  X  ( u  )  . 

Li  r 

In  considering  the  condition-method  formulation  F  (x  )=0.  which  is  independent 
of  any  potential  parameters,  we  have 

T. L,  r.  rj^/ 

F(x)=F(x{\i))  =  0, 

aF^/au®  =  (aF*7ax‘')(ax‘'/au®)  =  o  . 

The  last  equality  corresponds  to  the  relationship  between  the  condition  method 
and  the  parametric  method  for  the  same  adjustment  problem,  transcribed  in  matrix 
notation  as  BA  0 .  Here  this  equality  has  the  form 


With  regard  to  the  present  model  (67b)  and  its  geometrical  treatment,  we  have 
seeti  in  Chapter  4  that 


= 

P 


A  p 


N 


P  a 


R  R  ’ 

Thus,  the  above  equality  is  readily  verified: 


n'-a"  . 
r  a 


A  +  A  = 
a  a 


0 


6.2  Solution  for  a  Linear  Model 


In  a  linear  model,  by  definition  one  has 


A  ?  constant  /  0  ; 
a 


0 


0  . 


^  o  =  0  , 

a(3r 


The  equations  from  Chapter  4  presented  in  the  preceding  paragraph  are  restated 
for  the  insight  they  offer  in  the  linear  context: 

A®C>^  JfP*  fri^l  fA^l  ^ 

■  or  -  lA^l  . 


-A^R^  . 
A  K 


"k  '  4  • 


It  is  apparent  that  the  gradient  set  is  constant  and  the  model  surface  is  a 
hyperplane.  From  the  construction  of  the  Q-surface  It  then  follows  that  the 
latter  is  a  hyperplane  parallel  to  the  model  surface.  Further  equations  in 
Chapter  4.  such  as  (18a).  (19a).  ...  now  lead  to 


A  5  0, 
KX 


8  50,... 

Kxm 


and  equations  (18b),  (19b),  ...  lead  to 


N  5  0  . 
KX 


KXW 


0 . 


From  (41)  and  the  equations  below  it,  together  with  (39)  or  (39').  we  have 


%  -  V ^ 

SiiKc  (42)  and  the  equations  below  it  show  that  dx^dx^sO  for  any  dx^,  it 


follows  that  50 

Or 


Accordingly,  the  relations  preceding  (42)  yield 


A-  3  0 
i9r 


40 


with  similar  identities  holding  true  also  for  higher-order  partial  derivatives. 
From  equation  (43).  together  with  (38)  and  (34).  we  can  now  write 

.  a  .«a  .  a  .,L .  m 

Ax  =  Ax  =  <5x  -  g  hr  Q„,  N  <5x 

s  KL  m 

where  the  first  equality  is  consistent  with  the  fact  that  the  model  surface  and 
the  Q  surface  form  two  parallel  hyperplanes:  furthermore,  from  (44)  or  (45a, b) 
we  can  write 

.a  ,  a  r,®*  ^ 

Au  =  du  =  R-Ax 
P 

r 

The  remaining  subset  of  Ax  then  follows  from  (46)  or  (46')  as 


7.  CONCLUSION 


The  non-iterative  approach  undertaken  in  the  present  study  is  built  on  the 
geometrical  concept  of  the  functional  form  of  the  model  surface,  coupled  with 
the  new  concept  of  the  Q-surface.  The  geometrical  restriction  to  the  model 
surface  corresponds  to  the  restriction  making  the  observables  consistent  with  a 
given  adjustment  model.  The  fact  that,  in  general,  the  observational  point  Q 
does  not  lie  in  the  model  surface  corresponds  to  the  realization  that  actual 
observations  will  not  satisfy  the  adjustment  model  exactly.  From  the  geometric 
standpoint,  the  least-squares  principle  implies  that  the  desired  point  P  is 
the  foot-point  of  a  line  (straight  in  the  flat  observational  space)  dropped 
orthogonally  from  Q  onto  the  model  surface.  The  advantage  of  the  functional 
form  of  the  model  surface  stems  from  the  property  (Illustrated  in  three 
dimensions)  that  if  a  surface  is  expressed  by  N=constant.  then  the  gradient  of  N 
is  a  vector  orthogonal  to  this  surface.  It  is  then  possible  to  choose  a  point 
in  the  model  surface,  sufficiently  close  to  the  desired  least-squares  point  P. 
and  to  determine  the  amount  by  which  a  line  orthogonal  to  the  surface  at  this 
point  misses  the  "target”  Q. 

The  Q-surface  plays  a  crucial  role  in  the  determination  of  this  miss,  which 
is  expressed  in  the  Q-surface  coordinate  differences.  Such  differences  are 
ultimately  transformed  into  the  differences  of  model-surface  coordinates;  in 
adjustment  terminology,  the  latter  correspond  to  the  parametric  corrections.  An 
intermediate  transformation  features  another  set  of  model-surface  coordinates, 
linked  to  the  Q-surface  via  a  set  of  gradient  vectors.  All  transformations 
proceed  with  the  aid  of  the  appropriate  Taylor-series  expansions  containing  the 
first-  and  higher-order  partial  derivatives  of  the  observables  with  respect  to 
the  parameters  In  the  (nonlinear)  adjustment  model. 

The  outline  of  the  geometric  approach  to  a  one-step  nonlinear  adjustment  is 
presented  in  Chapter  3.  and  the  development  leading  to  an  algorithmic  resolution 
is  presented  in  Chapter  4.  The  centerpiece  of  this  approach  is  the  Taylor- 
series  equation  (43).  Its  first  term  contains  the  first-  and  second-order 
partial  derivatives  mentioned  above,  and  it  is  derived  in  detail  (see,  e.g., 
(!quatlons  39'  and  41):  the  second  term  contains  the  first-,  second-,  and  third- 
order  partial  derivatives,  and  it  Is  presented  without  a  derivation  (see,  e.g., 
equations  40  and  42);  the  third  term  contains  the  first-  through  fourth-order 
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partial  derivatives,  and  although  It  Is  not  presented  here.  It  Is  used  In  the 
numerical  example  of  Appendix  A;  higher-order  terms  could  be  developed,  but, 
similar  to  the  third  term,  would  be  useful  only  in  some  special  situations.  The 
varlance-covar iance  matrices  after  adjustment  (for  adjusted  parameters,  adjusted 
observations,  and  residuals),  and  other  adjustment  quantities  (residuals  and  the 
least-squares  quadratic  form)  are  presented  in  Chapter  5.  The  aspect  of 
reducing  a  nonlinear  adjustment  problem  with  nonlinear  constraints  to  a 
nonlinear  problem  with  fewer  parameters  is  presented  in  Appendix  B. 

In  closing,  we  summarize  the  main  features  of  the  non-iterative  approach 
dcv(>loj)ed  in  this  study,  and  compare  them  with  some  of  the  features  of  the 
"geometrical  method"  and  the  "extended  geometrical  method"  of  [Blaha,  1989a], 
which  are  iterative  in  nature: 

1)  The  present  approach  is  based  on  a  new  and  independent  concept  featuring 
the  Q-surface.  An  initial  point  is  not  displaced  by  iterations,  but  the 
solution  is  reached  in  one  step;  the  accuracy  of  th('  solution  depends  on  the 
number  of  terms  included  in  the  pertinent  Taylor-ser ies  expansions. 

2)  Although  the  algorithm  developed  herein  is  lengthier  than  the  algorithm 
for  the  geometrical  and  extended  geometrical  methods  of  [Blaha,  1989a].  the 
partial  derivatives  of  the  observables  with  respect  to  the  parameters  are 
evaluated  at  only  one  point  (the  initial  point  mentioned  above).  That  is  to 
say.  they  need  not  be  re-evaluated  at  new  points  as  is  the  case  with  iterative 
methods;  such  a  task  could  require  a  significant  effort,  depending  on  the  method 
and  accuracy  of  evaluation  (numerical  differentiation,  estimation,  etc.). 

3)  If  the  fourth-order  partial  derivatives  are  significant,  they  can  be 
incorporated  Into  the  derivations  and  calculations  upon  proceeding  with  the 
first  three  terms  in  equation  (43).  This  is  beyond  the  capability  of  the 
extended  geometrical  method  of  [Blaha,  19893],  which  can  utilize  up  to  tlie 
third-order  partial  derivatives,  and  far  beyond  the  capability  of  the 
geometrical  method  (ibid.),  which  can  only  utilize  the  flrst-and  second-order 
partial  derivatives.  (The  well  known  and  widely  used  linearized  approach,  where 
only  the  first-order  partial  derivatives  are  utilized,  may  be  severely  hampered 
by  the  nonlinearity  of  the  adjustment  model;  in  the  case  of  strong  nonlinearity 
coupled  with  inaccurate  Initial  parametric  values  and/or  large  observational 
errors,  the  solution  sought  in  an  iterative  piocess  may  diverge.) 
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4)  The  potential  advantagt;  of  the  non-iterative  approach  as  discussed  in  the 
preceding  item  Is  heightened  upon  considering  the  first  iteration  in  the 
extended  geometrical  method  of  [Blaha,  1989aJ,  where  only  the  first-  and  second- 
order  partial  derivatives  are  utilized;  it  is  further  heightened  with  respect  to 
the  first  iteration  in  the  geometrical  method  (ibid.),  where  only  the  first- 
order  partial  derivatives  are  utilized  (similar  in  this  respect  to  the 
linearized  approach).  In  some  cases,  the  first  iteration  could  actually  worsen 
the  Initial  parametric  values,  and  could  even  lead  the  iterative  process  astray. 

5)  In  following  the  theory  developed  herein,  it  is  possible  (albeit  tedious) 
to  extend  the  non-iterative  algorithm  so  that  it  could  eventually  utilize  even 
the  fifth-order  partial  derivatives,  which  might  be  of  Interest  in  some  special 
and  difficult  cases;  this  would  correspond  to  proceeding  with  the  first  four 
terms  in  equation  (43). 
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APPENDIX  A 


NUMERICAL  EXAMPLE 


The  present  numerical  example  Is  precisely  that  of  Appendix  A  in  [Blaha, 
1989a].  The  nonlinear  parametric  adjustment  model  is  represented  by  six 
observables,  each  expressed  by  a  third-order  polynomial  in  four  variables. 

Thus,  the  number  of  observations  is  six.  n=6.  and  the  number  of  parameters  is 
four,  u  =  4.  In  terras  of  the  four  parameters  u**.  a=l,2.3,4.  a  general  observable 
denoted  x  is  expressed  by 

a  ,  a  $  a  ^  Y 

x  =  c  u  +  k.u  u  +  m-  u  u  u  , 
a  a/S  afir 

where  a  ranges  from  1  to  4 .  fi  ranges  from  a  to  4,  and  y  ranges  from  /3  to  4. 

The  numerical  values  of  the  four  c-coefficients  (linear),  the  ten  k- 
coefficlents  (quadratic),  and  the  twenty  m-coef f Iclents  (cubic)  for  all  six 
observables  are: 
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.  1  . 

-0. 

.9, 

+  1  . 

8. 

-1 , 

.4. 

-0. 

9. 

+  2. 

1  . 

+  1 

.6, 

+  1 

.0, 

-2 

.  1 . 

+  1 

.  1 ; 

-1 

.7, 

+  2 

.2, 

+  1 

.6, 

-1 

.6, 

+  1 

.6. 

-1 , 

2, 

+  2, 

.  1 . 

+  2 

.3, 

-2, 

.  1 . 

+  1 

.2; 

-2 

.2. 

-0 

.8. 

+  1 

0. 

-1  . 

2, 

+0. 

9. 

+  1 

3. 

-0. 

8, 

+2. 

1 . 

+  1 , 

9. 

+  1 . 

2. 

+  0 

.9, 

-1 

•  5 , 

-1 

.6. 

-0 

.8. 

-1 

.0, 

+0 

.9. 

-0 

.9, 

+  2 

.  1, 

tO 

.8, 

-1 

.3. 

+0. 

.9. 

+  2. 

1  . 

-2 

.  1  . 

+  1 . 

8: 

+  1  . 

2, 

+  1  . 

7. 

+  2. 

2, 

-1  . 

5. 

+  1  . 

9, 

+  1  . 

3. 

-2. 

2. 

+  1 

.7, 

<■1 

.0, 

+  2 

.3: 

+  1 

.  1  . 

-1 

.2, 

-0 

.9, 

+  1 

3, 

♦^1 

.5, 

+  1 

.0. 

-0 

.8, 

-2 

.2. 

+0 

.8, 

+  1 

.7  . 

-1 

.2, 

-0 

,8. 

+  ]  , 

.  1 . 

-2, 

0. 

+  1 

.7. 

-0. 

9. 

+0. 

8. 

+  1  . 

3. 

-1  . 

8. 

-1 . 

1  . 
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+  1.7. 

+  1.2, 

+  1.7, 

-1.6; 

+  2.1, 

+  1.6, 

-1.1. 

+  1.9, 

+  1.6, 

+1.7,  +2,0, 

+  1.2 

+  1.6. 

-1.7; 

+  1.3, 

+2.1. 

-1.0. 

-0.8, 

+  1.2. 

-2.1. 

+2.2, 

-1.6, 

-1.8,  -0.8, 

H  .9 

+  1.2. 

-0.9, 

+  1.2, 

+  2.0, 

-1.7. 

+  1.1. 

+  1.0. 

+0.8. 

-2.2. 

The  positive-definite  variance-covariance  matrix  of  observations. 


corr(;sponding  to 

the  associated  metric  tensor 

rs 

g  .as 

chosen  as 

^+1.3266  +0.5254 

+0.3669 

+0.1572 

-0.2667 

-0.0549 

+1.3457 

-0.0281 

+0.5924 

-0.1448 

+0.5521 

fg  ]  = 

+1.4143 

+0.5144 

+0.3600 

-0.5204 

■ 

+1.4629 

-0.5604 

+0.4877 

+1.3930 

-0.1034 

+1.3500 

where  the  dots  denote  symmetric  elements.  The  positive-definite  weight  matrix 

of  observations  corresponds  to  the  metric  tensor  g  .  which  can  be  fovind  from 
ks  k 

g  g  -S  .  However,  this  tensor  is  not  needed  here, 
sr  r 

The  model -surface  coordinates  of  the  least-squares  point  P  are  c'^oscn  as 
u“  -  +1  ,4 .  -1.2,  +1.3.  +1.2  . 

where  the  notation  u*  could  have  been  used  Instead  of  u®.  The  observables  ”x" 

”  r 

for  this  point  are  denoted  x  ,  r  =  ],2 . 6;  computed  from  the  model  (third- 

order  polynomial),  they  are 

x*^  =  +4.7956,  -4.7901,  +33.4396,  +0.3331,  +22.4658,  -11.6899  . 

The  observational  point  Q  Is  generated  from  P  by  stipulating  the  values  of 
<5x'  .  where  <5x  is  the  vector  PQ  (see  Fig.  1).  Since  the  latter  is  orthogonal  to 

-  p 

tlie  model  surface,  the  (final)  residuals  correspond  to  -<5x  .  The  contravariant 

*  r  r  “ 

vector  <5x  is  generated  by  means  of  the  design  tensor  A  evaluated  at  P  as 

a 

follows  First,  the  metliod  of  Chapter  4  in  (Blaha.  1990]  gives  the  matrix 

FsfF'  F"]  grouping  the  contravariant  components  of  orthonormal  vectors  spanning 

r  r  r 

the  observational  space;  S.J  .  u  :  the  u  vectors  fi.  j,  ...  span 

the  model  surface,  and  the  n-u  vectors  v.  ...  are  orthogonal  to  it.  The 
r  r 

(  fimponents  of  C  ,  j  .  .  .  .  are  grouped  ante  the  columns  of  F'.  while  the 

r 

components  of  i/  ,  ...  are  grouped  into  the  columns  of  F". 
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Next,  llie  compoiuMits  of  <5x^  .  cori'ospondinK  to  obsorvat  i  oiia  J  errors  wliiih 
must  bo  recovered  by  *  he  adjustment,  are  generated  from  v*  ,  ...  as 

j-r  r  r 

<5x  =  au  f  h  71  ‘  ...  . 

IlcM'o  we  fiHve  n  ii  2 ,  and  choose  a  =  10.  b  -  35 ;  F"  (not  listed)  then  yields 

Sk  t3.21‘J3.  +22.5470.  10.5802.  -1.2758.  +2*.  2738.  +14.0805  . 

We  will  use  two  cases  of  generated  errors:  case  2  aj  above,  and  case  1  where 

a-1.0,  b=3.5,  and  wnere,  accordingly,  all  the  values  in  <5x'  are  scaled  down 

-2 

by  10.  The  least-squares  quadratic  form  corresponds  to  <5s  .  where  5s  is  the 
magnitude  of  the  vector  PQ.  In  general,  we  have 

2  --r  .-s  .  r.-  ?  ,2 

<5s  --  <5x  K  <5x  =  <5x  <5x  =  a  +  b  +  .  .  .  . 

rs  r 

In  the  present  example,  it  then  follows: 


case  1  . . . 

.  .  (5s ^  = 

2  2 

1  +  3.5  - 

13.. -'5  . 

( A  .  1  a  ) 

case  2  . . 

.  .  .  <5-s2  = 

2  2 

10  +  35  = 

1325  . 

(A. lb) 

These  values  must  be  recovered  as  the  weighted  sum  of  squares  of  residue  Is. 

"  r  “  r 

The  simulated  observations  correspond  to  x  +  <5x  .  determining  the  point  Q. 
However,  the  construction  implies  that  the  point  P.  obtained  by  the  orthogonal 
projection  of  Q  onto  the  model  surface,  must  be  the  same  for  both  cases  1  and  2, 

-  f 

as  wel  as  for  any  other  set  Sx  generated  as  above.  At  a  general  point  of  the 

model  surface,  the  model -sui face  metric  tensor  a  the  "necessary"  associated 

aj3 

rs  rs 

metric  tensor  g'  '  .  and  the  "ccmplementary"  tensor  g"  can  be  found  from 


A^'g  . 

a  rs  0 


(A. 2) 


,rs  ,r  a0  .s 
^  ■ 


(A, 3) 


rs  ,  rs 

g  ~  g 


(A. 4) 


where  the  model -surface  associated  metric  tensoi  a  is  such  that 


a/3  .a  r  I  r  1  1 

a  a.  =  5  .  or  (a  ]  =  [a„  J 

PS  r  s  PR 

If,  in  matrix  notation,  we  identify  (g  J=g*.  [g  '  ]  "g  ■  (g'  J=g'.  fg”  ]  =  g" 

fA  )=A,  (a  ]sa*.  and  [a  ^]=a,  where  we  have  attrltnited  to  (purely) 

<x  ctp 

covariant  tensors,  we  nave 


a 


a 


-1 


(A.Ca.b) 


* 


T 

A  K*  A 


(a*) 


g  '  A  a  A 


L'  e  L'’ 


(A  .  6a  ,  Li) 


If)  a  linear  adjustment  model,  a*  in  (A. 5a)  is  the  well-known  matrix  of 
normal  equations,  a=(a*)  ^  is  the  variance-covariance  mat’ ix  of  adjusted 
parameters,  g'  in  (A. 6a)  is  the  (singular)  variance-covariance  matrix  of 
adjusted  observations,  and  g'  is  the  (singular)  variance-covariance  matrix  of 

p 

residuals.  In  terms  of  geometry,  such  a  model  entails  constant  values  of  A 

q-A  P  <j  P3  pg 

and  thus  also  of  a  .,  a  .  g'  '  ,  and  g"  (in  addition  to  constant  g  and  g  ). 

ap  rs 

Ifi  a  nonlinear  model,  tliese  tensors  are  variable,  and  must  be  evaluated  at  P 
111  order  to  correspond  to  adjusted  quantities.  As  is  shown  in  Chapter  3  of 
{Clatia  1990J,  alternate  formulas  for  g'  and  g”  are 

g'  -  F'  ,  g"  -  F''F”^  ,  (A.7a.b) 

wiiere,  for  the  same  reason,  fF'  F"]=F  must  again  be  evaluated  at  P. 

All  of  Die  formulas  and  values  listed  co  far  have  served  merely  for 
verification  of  the  numerical  results  In  the  present  example;  the  actual  results 
hav  ’  been  computed  according  to  cnapter  5.  Thus,  when  applied  at  P,  the 

—  p 

formulas  (65)  or  (65')  have  recovered  <5x  for  cases  1  and  2,  and  equation  (66) 

r  s 

has  recovered  (A.la.b).  Similarly,  the  numerical  outcome  for  [g"  ]  computed  by 

(59)  or  (60a-c)  agrees  with  that  obtained  from  (A. 6b)  and  (A. 7b),  and  the 

p  g 

numcrica]  outcome  for  (g'  j  oompuled  by  (62)  agrees  with  that  obtained  from 
(A. 6a)  and  (A. 7a).  Finally,  [a  }  as  computed  by  (64)  agrees  witfi  (A. 5b),  while 
(a^^J  lifts  been  by-pa.ssed  (similar  in  this  respect  to  [g^^]).  Although  the  above 
matrices  are  associated  wltti  the  1  east -squares  point  P,  the  notation  Is  used 
wit(iout  an  overbar. 


The  variaiHC  covai  >ancc  matrix  for-  the  adjusted  parami.-ters  and  the 
variance -covariance  matrix  for  the  adjusted  observations  are  listed  as  follows: 


ft! 


«iS  , 


r+0. 009205 

i 

I 


+0.010983 

+0.022411 


-0 . 002919 
-0.010135 
+0.CO6404 


+0.005380 

+0.016375 

0.009214 


4  8 


+0.016113 


-0.5035  -*0.0501  -0.1940  | 

0.1372  -0.1325  *0.1852  | 

+0.0285  +0.8233  -0.4773 

+0.6167  -0.0587  +0.3701  ^ 

+0.6775  -0.2978  ^ 

+1.1668| 

There  is  no  need  to  list  the  variance -covariance  matrix  for  the  residuals.  Tha  , 

r*s  rs  rs  rs 

matrix  is  simply  [g"  ]  =  [g  ]  Ig'  '  1.  where  [g  ]  has  been  presented 

numerically  at  the  outset. 

As  in  Appendix  A  of  [Blaha,  1989a],  the  initial  set  of  parameters  is 
represented  by 

u®  =  +1.399,  -1.201,  +1.299,  +1.099  , 
o 

describing  the  initial  point  P  by  model-surface  coordinates.  The  observational - 

p 

space  coordinates  x  of  this  point  are  determined  from  the  model ,  as  are  the 

r  r  r 

sets  of  partial  derivatives  A  ,  fi  and  0  „  (the  last  set  is  constant  here, 

a  a/3  a^r 

and  any  further  set  is  zero).  The  above  coordinate  set  u^  will  be  utilized  for 
both  cases  1  and  2.  Accordingly,  In  the  system  (u**)  the  differences  between  P 
and  P  are  0.001  units  in  each  of  the  four  coordinates. 

We  already  saw  that  in  the  method  dubbed  "standard"  in  Appendix  A  of 
[Blaha,  1989a],  the  successive  iterations  diverged.  Below  we  recapitulate  the 
pertinent  results  only  for  the  methods  dubbed  "geometrical"  and  "extended 
geometrical"  (ibid.).  In  the  former,  we  present  the  worst-case  errors  in  the 
first  three  iterations,  while  in  the  latter,  we  present  such  errors  in  the  first 
two  iterations;  errors  In  subsequent  Iterations  are  completely  negligible. 
Immediately  following  the  reviewed  values,  we  list  errors  resulting  from  the 
non  iterative  approach  developed  in  the  present  study.  Three  results  will  be 
listed  for  cases  1  and  2: 

(a)  The  worst-case  error  in  u**  stemming  from  the  utilization  of  only  the  first 
term  on  the  right  hand  side  of  (43).  linked  to  the  development  of  equation  (41); 
this  error  is  referred  to  as  "1-term". 

(b)  The  worst-case  error  in  u*  stemming  from  the  utilization  of  the  first  two 
terms  on  the  right-hand  side  of  (43),  linked  to  equationr  (41)  and  (42),  of 


which  the  latter  has  been  listed  without  derivation;  this  error  is  referred  to 
as  "2  term" . 

(c)  The  woist  case  error  in  u*  stomroing  from  the  utilization  of  the  first  three 
terms  on  the  right-hand  side  of  (43),  of  which  the  third  term  has  not  been 
presented  at  ail;  this  error  is  referred  to  as  "3-term". 


According  to  the  above  description,  the  results  are  presented  as 


case  1 

geom. : 

4 

-2x10 

.  +1x10 

o 

! 

o 

X 

1 

ext .  geom . : 

-3x lO'® 

,  -Ixio" 

10 

non- i ter . : 

1  term ; 

+  2x10 

"8 

2-term:  +4x10 

3-terin:  +1x10 

case  2  , 

geom . : 

-2x 10  ^ 

ixio” 

5  -8 

,  +1x10  ; 

ext .  geom . : 

-lxlo“^ 

.  -2xl0“ 

10 

non  -iter . : 

1  -  term : 

+3xl0'^  . 

“8 

2-tcrm:  -4x10 

3 -term;  -5x10 

An  additional  variation  of  the  present  example  has  also  been  resolved, 

where 

u“  =  +1 .39,  -1.21,  +1 .29,  +1.09  . 
o 

The  differences  between  P  and  P  in  the  system  (u**)  are  now  0.01  units  in  each 
of  the  four  coordinates,  l.e.,  ten  times  larger  than  previously.  This  set  u* 
is  utilized  ofily  in  conjunction  with  the  observational  errors  of  case  1.  The 
results  can  be  imagined  written  to  the  right  of  the  above  case  1:  this  would 
.illow  one  to  visualize  the  growth  of  the  worst-case  errors  In  u**  due  to  a  ten¬ 
fold  increase  in  observational  errors  (by  looking  downwards,  l.e.,  examining  the 
above  case  2),  or  due  to  a  ten-fold  increase  in  "errors"  in  the  initial  values 
(by  looking  to  the  right).  The  results  stemming  from  this  last  variation  are 

-  3  -5  -7 

case  1  ...  geom..  3xl0  '  .  »lxl0  ,  +lx]0  ; 

ext.  geom.:  2x10  '  ,  +4x10  ; 

non-iter.:  1-term:  +2xlo  ^  .  2  term:  +3x10  ^  ,  3-term:  +2x10  ^  . 


In  all  three  methods  ("geom.",  "ext.  geom.",  and  "non  iter."),  the  increase 

In  errors  in  the  Initial  values  appears  to  have  much  moi  serious 

o 

repercussions  than  the  Increase  in  observational  errors.  In  fact,  the  results 
in  "ext.  geom."  and  "non-iter."  seem  to  be  degraded  only  slightly  by  the 
increase  in  observational  errors.  Iterations  in  "geom."  and,  especially,  "ext. 
geom."  are  seen  to  decrease  the  worst-case  errors  in  u*  more  dramatically  than 
the  addition  of  terms  in  "non-iter.".  In  the  latter  case,  the  order  of 
magnitude  of  such  errors  is  seen  to  decrease  fairly  regularly  by  an  Inclusion  of 
additional  terms,  consistent  with  the  characteristics  of  an  expansion  in  the 
Taylor  series.  However,  the  above  results  arc  skewed  in  disfavor  of  the  non¬ 
iterative  method,  in  that  only  a  third-order  polynomial  serves  as  an  adjustment 

I' 

model.  In  such  a  case,  the  third-order  partial  derivatives  grouped  into 
are  constant,  and  the  extended  geometrical  method,  which  cannot  utilize  the 
fourth-  and  higher-order  partial  derivatives,  is  spared  any  harmful  effects 
these  derivatives  could  entail.  By  comparison,  the  third  term  (not  presented) 
in  the  non-iterative  method  operates  with  the  fourth  order  partial  derivatives, 
but  the  latter  arc  identically  zero  here. 

We  can  conclude  the  discussion  by  stating  that  if  the  fourth-order  partial 
derivatives  of  the  observables  with  respect  to  the  parameters  were  significant, 
the  quality  of  the  results  produced  by  the  non-iterative  method  would  likely 
improve  vis-a-vis  the  extended  geometrical  method.  The  Improvement  would  be 
more  dramatic  vis-a-vis  the  geometrical  method,  which  utilizes  only  up  to  the 
secono  order  partial  derivatives.  Such  an  improvement  would  manifest  itself  in 
yet  another  respect.  This  respect  comes  to  light  upon  considering  the  very 
first  iteration  in  the  extended  geometrical  method,  in  which  only  the  first-  and 
second  o’'der  partial  derivatives  are  utilized.  Accordingly,  the  extended 
geometrical  method  presents  potential  danger  that  the  first  Iteration  could 
actually  result  in  worsening  of  the  initial  values  u*  in  some  special  cases, 
and  could  even  lead  the  iterative  process  astray.  With  regard  to  the 
geometrical  method,  such  potential  danger  would  be  heightened  because  this 
method  treats  tlie  first  Iteration  in  a  linearized  fashion,  l.e.,  utilizes  only 
the  first-order  partial  derivatives. 


APPENDIX  B 


PARAMETRIC  ELIMINATION  DUE  TO  NONLINEAR  CONSTRAINTS 

B . 1  Initial  Parametric  Values 

In  standard  adjustment  notation,  a  set  of  s  constraints  "G"  among  u 
parameteis  is  symbolized  by 

G(X'’)  -  0  , 

where  represents  the  set  (column-vector)  of  adjusted  parameters.  The  usual 
approach  consists  in  expanding  these  constraints  In  the  Taylor  series  using  an 
initial  set  of  parametric  values,  symbolized  by  X*^.  It  then  follows  that 

G  ( X^ )  =  W  C  X  +  .  .  .  =  0  . 

C 

where 

W  =  G(X°)  .  C  =  OG/aX)  . 
c  o 

and  where  X=X^-X°  is  a  set  of  parametric  corrections.  The  subscript  "o" 

indicates  that  the  matrix  C  is  evaluated  using  the  elements  of  X° .  This  matrix 

has  the  dimensions  sxu,  while  the  column  vectors  W  and  X  contain  s  and  u 

c 

elements,  respectively  In  a  standard  (linearized)  approach,  the  terms 
symbolized  above  by  dots  are  omitted. 

Here,  as  well  as  in  standard  adjustment  theory,  the  constraints  are 
considered  independent,  in  the  sense  that  the  constraint  matrix  C  has  the  full 
row  ranl<  s.  However,  the  standard  theory  proceeds  in  general  with  By 

contrast,  the  current  nonlinear  development  will  benefit  from  such  initial 
values  X°  for  which  it  holds  true  that 

G(X")  =  0  .  (B.l) 

In  this  case,  only  u-s  of  the  u  values  in  x”  can  be  chosen  independently. 
Althougli  a  set  X°  for  which  G(X^)/0  would  also  he  acceptable,  the  computations 
below,  which  lead  to  G(X°)“C  would  eventually  have  to  be  performed  as  well,  and 
the  resulting  formulas  would  be  more  cumbersome  and  less  tractable. 
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We  present  an  example.  In  matrix  notation,  of  a  simple  iterative  algorithm 

leading  to  G{''°)=0.  First,  the  matrix  C  is  partitioned  into  [C^  C^]  .  where  the- 

submatrices  and  have  the  dimensions  sx(u-s)  and  sxs.  respectively.  The 

submatrix  can  be  considered  nonsingular  without  any  loss  of  generality, 

since,  if  needed,  tlie  parameters  could  always  be  rearranged  beforehand  for  this 

condition  to  be  satisfied.  The  vector  is  similarly  partitioned  into  the 

subsets  denoted  temporarily  Xj  and  X^.  which  contain  u-s  anil  s  elements, 

respectively.  The  elements  of  X^  are  chosen  and  held  fixed.  The  remaining 

elements,  grouped  in  X  are  subject  to  change.  We  symbolize  their  initial 
0  ^ 

choice  by  X  to  accommodate  the  iterative  indices  below.  The  first  iteration 

1  1 

yields  corrections  grouped  in  AX^,  resulting  in  i.n  improved  vector  X^.  After  an 
i-th  iteration,  the  improved  vector  becomes 


X 


1 

2 


+  AX 


1 

2 


If  the  corrections  become  negligible  and  the  iterative  process  is  terminated 
after  n  iterations,  the  values  in  X^  are  adopted  as  the  elements  of  X^.  In 
joining  this  vector  to  the  independently  chosen  vector  X^ ,  one  obtains  the 
desired  vector  X^’ . 


The  model  before  the  1-th  Iteratin',  can  be  presented  as  follows; 

G(X°)  =  GlXj.X^)  =  GiXj.X^'^)  +  (Xg  -  X^'S  +  ...  =  0  . 

resulting  in  the  1-th  solution; 

"a  -  4'’  •  '  '■■IX, -4'’'  • 

Upon  neglecting  the  higher-order  terms  represented  by  the  dots,  one  writes 


.  4  ^  x‘  ‘ 


(C^  ^  GlXj.Xg'^) 


If  the  values  in  AX^  are  not  deemed  negligible,  the  matrix  C  is  updated  upon 
evaluating  0G/8X  with  the  new  values  X  and  is  denoted  C  .  Similarly, 
one  forms  an  updated  vector  GIX^.X^)  and  proceeds  to  the  iteration  i+1.  After 
till-  final  x”  has  been  computed  by  this  or  a  different  method,  an  actual 


nonlinear  least  sijuares  adjustment  with  constraints  can  take  place. 
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B . 2  Model  Subsurface 

As  wo  have  seen  in  Chapter  2.  the  nonlinear  parametric  adjustment  model  can 
be  written  as 

'',a.  *r.o 

X  (u  )  =  X  u  )  ^  A  Au  +  (1/2)  11  .  Au  Au  ♦ 
o  a  a/3 

(1/6)  •t’'-  Au”  Au^  Au^  +  ...  .  (B.2a) 

o|3r 

where 

Au”  =  u”  -  u”  .  (B.2b) 

o 

with  r  =  1.2 . n,  a  =  1.2 . u;  and  where 

,r  .r  /^2r_a-0. 

A  =  Ox  /0u  )  .  Q  -  =  (9  X  /au  9u  ) 

a  o  a/3  o 

/^^r,.  a_/S^r. 

<t>  „  -  (9  x  /9u  9u  9u  )  . 

aBr  o 

All  lower  case  Greek  letters  vary  in  the  fashion  prescribed  above  for  o.  A 
similar  convention  applies  for  other  kinds  of  indices  as  well  (lower  case  Roman 
letters,  etc.).  The  subscript  "o"  indicates  the  evaluation  at  the  initial  point 
P  lying  in  the  model  surface. 

The  s  nonlinear  constraints  Joined  to  this  model  are  represented  by 

G^(u“)  =  0  ,  (B.3a) 

where  L=l,2 . s.  This  equation  can  be  regarded  as  the  functional  form  of  a 

surface,  generalized  to  higher  dimensions.  In  referring  to  (B.l).  one  also  has 

G^{u®)  =  0  ,  (B.3b) 

where  the  values  u”.  a=],2 . u.  are  known.  The  functional  form  (B.Ca) 

o 

restricts  the  final  least -squares  point,  whose  model -surface  coordinates  are  u  , 
to  a  certain  lower  dimensional  surface  embedded  in  the  model  surface.  From 
(B.3b)  it  follows  that  the  Initial  point  P  also  belongs  to  this  lower- 
dimensional  surface,  which  will  be  called  "model  subsurface".  This  is  apparent 
from  the  right-hand  sides  of  (B.3a.b),  which  contain  the  same  sets  of  constants 
( zeros  ) . 
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The  current  development  is  organized  along  the  foiJowing  lines.  First,  the 

coordinate  set  {u*),  a=].2 . u.  is  partitioned  into  {u^.u*^},  A=1.2 . ii~s, 

K^l,2 . s.  This  allows  (B.3a)  to  be  written  as 

u*^)  -  0  .  (R.4  ) 

representing  the  functional  form  of  the  model  subsurface.  Subject  to  the 
condition  stated  explicitly  in  the  sequel.  (B.4)  makes  it  possible  to  express 
the  last  s  coordinates  In  terms  of  the  first  u-s  coordinates; 

u=u(u).  (B.5) 

Kquation  (B.5)  is  the  Monge  form  of  the  model  subsurface  embedded  in  the  model 

surface.  wlu;rc  u^.  A=1.2 . u  s.  are  the  subsurface  coordinates  (independent 

variables).  The  substitution  of  (B.5)  into  (B.4)  yields 

g'^^(u'')  3  G^(u''.  u'^(u'')  )  =  0  .  (B.6) 

whicli  is  an  identity  in  the  model  subsurface.  Thus,  further  identities  follow; 

L,.  A  .  .2  L..  A-  0  .  -3  L,^  A.  0.  f 

dg  /  du  =  0,  d  g  /  du  du  =  0.  d  g  /dv  dv  dii  =  0 .  (B.7) 

K  A 

which  will  be  seen  to  lead  to  a  relation  for  0*^(0“).  In  this  manner,  the 
parameters  u  .  K=1.2 . s,  will  have  been  effectively  eliminated. 

Expressed  in  the  Taylor  series.  (B.3a)  reads 

G^^(u*)  =  Au*  +  (1/2)  Au^  Au^  +  (1/6)  T*'^  Au*  Au^  Au^  +  ...  =0  . 

a  «Pr 

wtiere  advantage  has  been  taken  of  (B.3b).  and  where 

=  (aG^/au“)  .  h''  -  (a^G^/du®au'^)  . 

a  o  a0  o 

=  (a^G^/au“au*^au’’)  . 

afir  O 

In  using  the  partition  of  (u*)  and  the  symmetry  of  partial  derivatives  in  the 
lower  Indices,  one  develops  this  equation  as 

L-  A  K.  ^L.A  a  ^  O  tr^  A  ^  A  ^ 

G  (u  .u  )  =  Au  C,.  Au  +  (1/2)  (H._  Au  Au  ^  2  Au  Au 
A  K  AU  AK 

Au*^  Au*^)  +  (1/6)  Au^  Au***  +  3  T^^^^  Au^  Au^  Au*^ 

^  3  Au^  Au*^  Au^  +  Au*^  Au^  Au^)  +  .  ,  =  0  .  (B.8a) 
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where 

Au=u-u.  Au=u-u.  (B.8b) 

o  o 

and  where  u^.  u*^,  the  mode] -surf ace  coordinates  of  P.  are  known  (see  the 
o  o 

relation  B.3b  and  the  statement  below  it).  Equation  (8.8a)  corresponds  to  the 
step  represented  by  (B.4).  In  evoking  (B.5),  we  next  formulate  the  Taylor 
series  for  u  : 


Au*^  -  A'^  Au^  +  (1/2)  Au^ 

A  All 

where 

.  ,  K  , ,  K , ^  A, 

A.  =  (3u/au  ) 

A  o 

The  partial  derivatives  are  again 


.  /A  a  a  ,  n  w 

Au  +  (1/6)  ®  Au  Au  Au  + 
AllV 


Q' =  (a  u  /au  9u  ) 

AQ  o 

^,K  _3  K,^  A^  0^  <V, 

®  AOV  = 


symmetric  in  the  lower  indices. 


(B.9) 


The  substitution  of  (B.9)  into  (B.8a)  in  view  of  the  step  (B.6)  yields 
after  a  convenient  rearrangement: 

,L  „L  .  ,K  >  .  A  _  ,„L 


L,  A, 
g  (u  ) 


(C“  .  C“  A-',  AU-  .  (1/2,  (C-  , 


2  A-*^ 

AK  n 


^  »KM  "'a  ®'aq4.  ^  ^ 


.OIJ^  a|K  ,M  .>pf'  *•»  A '  ^ 

^  ^  A  ”  ni?  ^Anv  ^  ’^AQK  ^  »P 


+  3  T^  A-*^  A-” 
AKM  O  f 


„I-  A  .*<  A  |M  .  .N  ,  ,  A  ,  0  ^  <P 
^  ^KMN  ""  A  ^  n  ^  " 


=  0 


This  identity  is  immediately  confirmed  at  P.  Upon  differentiating  it  in 

6  0  r 

succession  with  respect  to  u  ,  u  u  ,  etc.,  and  then  rearranging  the  free  and  the 
dummy  indices  while  taking  advantage  of  the  symmetries  in  partial  derivatives, 
in  accordance  with  (B.7)  it  follows  that 

a,‘-/3u''  .  c);  .  A'l;  ,  (c|;  •  b;;„  ,  A'*  .  (4 

•  "in  “'a)  ♦'aO*  •  "m  •  4* 


+  2  T^'  A'*^  +  T^  A’*^  A’^  +  2  (H^ 

AOK  f  AKM  n  V  '  QK 


A'”l 

o'  ^  AV 


f  (h‘'  1)-”  +  t'^  .  2  T^  A'” 

“km  “  nv  ^  ^OKM  ^  f 


_L  .,M  ,,N.  *K-  ^  0  ^  I'  „ 

^KMN  0  ''aJ  .  ...  =  0  , 


(B.IO) 
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( [i .  n  ) 


^2  L,.  A,  n  ^,K  ,,1.  „L  .,K 

a  g  /au  au  =  a'^  . 


+  { ID '  '^  t  f  .  .  .  1  }  All'*’  * 

'  K  AdW  '■  ‘  ' 


K 

A 


-  0 


a^g7au^au**au'*’  =  M  •••  )  ^  ••  -  =  o  ;  (b.12) 

the  expression  [  ...  ]  appearing  in  {B.11,12)  is  presented  explicitly  below  as 
the  expression  inside  the  brackets  in  (B.15). 


The  evaluation  of  (B. 10-12)  at  P  yields,  respectively, 
,  J 


A' 


O' 


A 

J 

AO 


”7;: 


"l  '"ao 


*.«  /uL  „L  .,M,  ..K- 

“aK  ^  0  ^  **0K  “km  ^  0  ^  A  ^ 


(B. 13) 


(  B  .  1  4  ) 


where 


<D'*  =  I)‘*  [H*'  0'*^  +  T*"  1  T*"  A'*^  *  T*'  A'*^  +  T*'  A'**'  A'''* 

AOV  L  ^  AK  Of  AOf  AOK  f  AfK  0  AKM  O  f 


*  «K« 

A•^ 

0  ’ 

0'*^  + 
Af 

(H*'  + 

'  fK 

H*'  A'”) 

"km  f' 

O'*" 

AO 

^  <»KM 

+  T*" 
OfK 

+  T*^ 
OKM 

A-”  . 

tL  A'” 

fKM  0 

"  ^KMN 

A.> 

'1)  A- 

^  • 

(B.15) 

h  S  -  \ 


The  last  equation  represents  the  condition  mentioned  below  (B.4).  In  matrix 

M  L 

notation,  this  condition  states  that  the  matrix  [d”]  Is  the  inverse  of  [C*"). 

L  L  ^  ^ 

which  in  turn  implies  that  the  matrix  [C  C  ]  of  dimensions  sxu  must  have 

A  K 

the  full  row  rank  s.  and,  therefore,  that  the  constraints  must  be  linearly 
independent.  In  the  affirmative,  an  eventual  rearranging  of  parameters  will 
ensure  that  the  matrix  [c|^]  is  regular  (i.e..  square  and  nonsingular).  This 
subject  has  already  been  discussed  in  Section  B.l,  and  has  led  to  (B.3b).  Upon 
substituting  (8.13-15)  and  higher-order  partial  derivatives  (not  listed)  into 
(B.9),  one  obtains  a  relationship  for  u  as  has  been  indicated  below  (B.7). 
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B . 3  Modified  Nonlinear  Model  with  a  Reduced  Number  of  Parameters 


Having  expressed  the  model  subsurface  in  (B.5)  via  (B.9),  we  can  reduce  the 
number  of  parameters  in  the  model  equation  (B.2a)  by  substituting  into  it  the 
expression  for  Au  from  (B.9).  Upon  the  realization  that 


.r  .  a  .r  ,  A  .r  .  K 

A  Au  =  A.  Au  +  A„  Au 

a  A  K 


„r  .  a  .  $  ^r  ^A^ft  ^  .A.K  *K.M 

0  a  Au  Au  Q,^  Au  Au  +  2  0.„  Au  Au  +  Au  Au 

a3  AQ  AK  KM 


-  (3  .  r 

®  .  An  Au  Au 

a^r 


.r  *A,0^4'  „.r  A^O^K 

'  Au  Au  Au  *  3  Au  Au  Au 

AOU’  AOK 

o.t.f  *AKM  .r  .K.M.N 

+  3  <!».„„  All  Au  Au  +  Au  Au  Au 

AKM  KMN 


this  substitution  yields 


c  =  x^  .  (A^  .  A'^)  AU  .  (1/2)  (n^^  .  A^  0 

♦  2n^  A'*^  +  n*^  A'^A'^)  Au^  Au^ 

AK  fi  KM  A  0 ' 

(1/6)  +  Aj^  ‘‘’'aqip  ^  ^AK  ^  ‘^KM  ^  A 

.  ,K  o  .  tK  .  ,M 

4-  3  (1)  A'  +  3  d)  A*  A* 

AOK  V  AKM  *  n  »P 

.  .X  *  .M  ^  ,N  ,  ,  A  ^  0  ^  4> 

^  ®KMN  ""  A  0  ^  4^^  ^  •••  • 


.K 


(B.16) 


,  K 


,  are  known  from  (B. 13-15).  etc.  The  symbols  x  and 


where  A '  _  O' 

A  Ad 
r 

x^  in  (B.16)  are  interpreted  as 

r  r.a.  r.A  K.A..  -r.A. 

X  =  X  ( u  )  =  X  ( u  ,  u  ( u  ) )  =  X  ( u  )  . 


r  r.a.  r.A  K.A.,  -r.  A, 

X  =  X  (u  )  =  X  (u  .  u  (u  ))  =  X  (u  )  . 

0  0  o  o  o 

where  is  known, 
o 

Equation  (B.16)  is  now  reformulated  to  read 

-r.  A,  -r.  A,  rr  ,  A  ,,  rr  »  A  ,  0 

X  (u  )  =  X  (u  )  +  A,  Au  +  (1/2)  n.„  Au  Au 
o  A  AU 


-r  A  0  V 
+  (1/6)  M,  Au  Au  Au  +  ...  , 
AQ4> 


(B.17) 


where  the  sets  of  implicit  partial  derivatives  at  P.  namely 
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'  (3x^7  3ii^) 


,  ^2  -r  .  -  A  ^  li , 

(3  X  7  3 11  3u  ) 


o 


/:i3-r,-  0^  ij'. 

®AQW  "  ^o 


follow  readily  from  (B.16): 

rP 


'A 


a''  a*'  a  .K 

^  ^  ^  A 


AO 


‘^Afl  "  '^k-  “'aO  '  2  “aK  '^’o  '  “km  ''’a  ^'o  • 

Afl  K  AO  AK  0  KM  A  0 


(D*^ 

AOI' 


0*^  +  A^  0'*^  +3  0*^  0'*^  +  3  o'^  A'*^  0'^ 

AOf  K  AOW  AK  O'?  KM  A  O'? 


o  as''  a.X  o  Asf  a.K  ..M  .r  A.X  a.M  ,,N 

♦  3<Ii  A  +3<1)  A  A*  +  0  A  A'  A 

AOK  <?  AKM  O  '?  KMN  A  O  '? . 


(B.  18a) 

(B. 18b) 

(B. 18c) 


Wtiereas  equation  (B.2a)  represents  a  nonlinear  model  in  the  parameters  u^ , 

a  =^1.2 . u,  equation  (B.17)  represents  a  nonlinear  model  in  the  parameters  u  . 

A-^1.2 . u-s.  In  geometrical  terms,  u*  are  the  model -surface  coordinates  of 

the  point  depicting  the  unconstrained  least-squares  solution,  and  are  the 
model -subsurface  coordinates  of  the  point  depicting  the  constrained  least- 
squares  solution.  Either  solution  can  be  carried  out  using  the  geometrical 
algorithm  described  in  the  body  of  the  present  study,  or  the  algorithm  from 
(lilaha,  1989a|.  However,  having  eliminated  the  s  parameters  u  ,  one  now  has 
a  smaller  system  to  resolve. 


In  conclusion,  in  the  parametric  adjustment  model  with  nonlinear 
constraints,  a  point  representing  the  least-squares  solution  Is  restricted  to 
lie  in  the  model  subsurface,  i.e..  a  surface  of  smaller  dimensions  than  the 
model  surface  in  which  it  is  embedded.  The  geometrical  approach  leads  to  the 
replacement  of  the  model  surface  by  the  model  subsurface,  and  to  the  treatment 
of  the  observational  point  Q  with  respect  to  the  new  surface  in  the  usual 
manner.  Accordingly,  the  constrained  least-squares  point  is  the  result  of  an 
orthogonal  projection  of  Q  onto  the  model  subsurface.  In  the  adjustment 
terminology,  this  approach  eliminates  s  of  th«!  original  u  parameters,  where  s  is 
also  the  number  of  constraints.  The  remaining  parameters  are  resolved  by  the 
method  of  the  nonlinear  parametric  least-squares  adjustment  without  constraints, 
where  all  the  arrays  must  be  properly  modified. 
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